


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations l. Thesis and Dissertation Collection, all items 


1988 


A study of second and third order models for 
the tracking subsystem of a radar guided missile 


Williams, John Walter. 


Monterey, California. Naval Postgraduate School 


http://hdl.handle.net/10945/23127 


Downloaded from NPS Archive: Calhoun 


Calhoun is the Naval Postgraduate School's public access digital repository for 


¿ К D U DLEY research materials and institutional publications created by the NPS community. 
| Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 
th 
KNOX appointed — and published — scholarly author. 


ІШ LIBRARY Dudley Knox Library / Naval Postgraduate School 


411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 








http://www.nps.edu/library 


“ 
1 
б 7 É 
Ы Ы C E 
. E 
Ш rm 
.4 АП 
. 
П 
on 
` 
... 
, 
П 
В О . E 
GE 
P 
О 
О 
1 0 E 
r 
E а "D ' 
О 
"E E 
О 
О СИМ 
П 
О . . 
LI 
a D 1 
Й 
TIRO te 
А б 
А р 
, 
б 
О 
О .а n . 
О 
б тах 
.. ПЕС .. 
|| ' 
n D Й 
Й 1 
r 
П 
. t» 
Й 
О е 1 
LEE . | . 
. . С Й ифа ' 
СИМ D 
^ P 
En] E 
Ú 
О . 
» P б А 
D r 
a , г ІШТІ! 
D . а 
te 
Й "P" 
g [ELT 
Й а 
. 
. ao б 
. ПА 
s [ELI 
» s П 
П E ПИО 
. а 1 . О 
. « bras 
е LI 1 
s 
. 
I 3 О С 
P 0 А 
П E 
» ‘ 
О D 
L] М 
ПИК] ' 
one LIBET О 
П ПЕСИ 
Li » Її 
ы E П А 
' П 
" а ри м до 
О 
ae | .  . 
М 
П » Tm 
б б 
L » 
Й П 
Й ` 
L ЛЬО О L] 
О П 
П P 1» 
n ' 
P 
Ы О 
4% toe 
1 
П Й 
б . 
E Й . 
П О 
' ПО П 
0 і 
4 
га 
D 
a 8 
P 
t 
» 
E 
, . 
[i ae 
П 
g 
* ) 
P 
.. " 
LAC) 
О 
О 
. 
4 . 
Й 
n б 
» 
. П 
Ы , y ET 


|. . 
A 
р 
E A 
В 
ПИО а 
. 35 
E AUD 
LET . LI 
A Е 
E TE 
` 
vs " 
а! LI 
“! ы О 
“ 
' 
а 
n Й 
+ LI 
Й 
J | 
ва 
. ~ 
5e ГЕ 
E 
E А ta 
. LI Ld 
"1 
* 
T n а 





nudam 7 U 
















































































































































"LL > =. г, 
з hi A A EAM LA А УУ Y dans y cia dota к ИЦ БА 
DN ON к T | 4» 3 за "A A k L tr Pire te DATAR ENT] КИЕ, ГИР ТЭР A ада, СИ ТАР 
n ә ж е, ; а, мера КУР Ру м ee OS АС ^re е 
а i E [1 4 2 а 14.14 t ААТ" аа a veri ОТЕ ады v Ae ym А 
t б ts ШИИ va E ье С РИС О И E О ТРОА У TOEN т C И el УДАА ty = МСД = = 
. r РА E в. "ann vw vt n раје“ 07 4л э! ment o Ра Ри ПС офа a 
D D " h ИСТО УТС rwn DDR MNT TT: и EM AY ARR и CODO! UG load > 
в a DE ‚у or ык э ғ” ПИ ДИН ПАДА а а УЋИ ИС ПУН Ferro pr rra XP E d Ln 
LN " A. TAN АМА ТАСТ ТУЗ DEDOS EU S DNUS TEMA пада и: эрге, ow 5 » n= 
n B . Си А | ғ. S T L А ^ DES K i Us een миа NUS rry peche МУ одаје у ча LRL POE roa к ies пи аи 
5 Ч oe * Пим а B xL и 4 <. Au [К ШЕТУ ҮҮ 4 Q n Я 4 
., ГЕН о d , [УРЕДИ ал КАЛТАТ ДУУ; MATES я PDT rh co riim "xr od НТ DE АЦ НАХ - V m opem ти 
А а “ә № d аА ТАЛ КИК КАЛД И. КЕ pis ы ТНҚ ең РОКИ эк А-аа ө. А Ала реінде е мачем оца 
мш СИ . о ЗЕ CD WD TO TR nnt AL A А [ro n и ao IH n P аа АК өр ae сафол аја 
В А . "P “<. af aa v" ETE MEER Ree) AAA a E S: кР УА ТҮТҮК ИП der tme Зе be ru E эләр чь тэу kaspaya 
' Г » РР 4 УТИ? ce ма! Ase а) пече Banca! КЫЗМАТЫ ЗИ, Y DUM Ту mU ПТ н 
. Y " Й О ОТЕ ЧУ РОСА: РУКА У РАСУ refers es pe DOE LEN еура Dato c > A ra 
0 ПШ `. П) à 'Aw e г ТЖЕ PP NOA t А.» Pa reo 42. ИДА) "FETTE" An de rers Ur MA оба чин d ар љето с ТО ЕД У dr 
Ç 6.) m TEM [E 4 о ага а re БЕЛО y PP OT LECCE. D ОНА A кои Мон пуни vp др потоа а rer erem 
aig 5 0 "A 1 1 Vi NM c (27 a UTA еги | ДЕЛ УК ОР АТИ МАГА а ју МИ И КА: MAN S dea Марини та dew Leda. mel. y А СТ 
d ' у ho 1.220.632 cuui y А m AT DULCE A CR ға S de Nn b аром Е Үе С D EE rd 
. ' ay .. Й Ж de 3 D" Б № я ae С 9 еа lal Ре ES Б eret О AT "t Aa M Мб PAM 2 s 
TTL У ЧИ СУ TA or А ғ 0,7 A a T Ma LETS у ЫН EE EM PLI я а P sim БЕ га” 
mn ова а a No z А л т р. ЖАПСА а CN AE WI LM YE win ей: Er oM үлүү ee МАН ^ Ве pee v у ow. За x 
.. А r с ПА AeA yo a A ҮҮ УГ AA A CETRO n A да ни SANTA ы», 80 
E Й E к r: а р СО" nu (METTRE! eG Hr O] ща И "n LINT T Phy topi Me nde y Ce Ou o e ы š 
О О ^ шуы ' м fe .а ДИУ AA AS Да 2. сета Е, TEC pre 04 е Dr AO т р Dos ress " ДИМ Uere do Mori 
Ë u * а si Š Aue “, 5 ыы У БУНА ИК АЧ ml oh e eri e а ута фуъл:т а ПУХ hiba ita ТЕЕ 
Los O к о В р " "T = 4 Xm КӨТ А лүүчү ЕМА DOE Му ы! КА че PARE ТУ) RR К TUUS + At 
B à " Џ " " ud 1 ка Й ЕРГЕ У | „1% у и Gard A T iv ма nw гута as Мот CP Ces LLL има E Alan Ma M 
E P "n" зе м дамо ИУ Fear "| митът чука? ка мк DOE d re “ir аи е кч 
, i М 4 ` E П уу es wi DP У YT y^ INTERN үзі РОС ТЕК ПА рен bs E Ld Lo Pep uris ар сев 
of @ 4 4 а Sy fae ГҮҮ, ñ y Pom ELA La хи кь» Aer Ч ТРТУ га" mm NG Ka gi hyna у» CIA DENM, m adi. уы за ма DU AAA apna 
А U " П 4 a ре "EON re IC ^A a ксле кү РИ И А ee ee ew ve la) mer ВЕРУ y [Am дер қ p ios абы. ee кант EI ы 
П “.. алыл Ф ча PUMA A E A TR L P LOL AT Pere i AM Пад Ll AAA, ARA ли 
Pr р E ' П 2... -ч... o ER AN WP TT WP ЛАЗ tariy > Поатие ПЕТУ КОРИ A У ДНЕС мр тиара ГА 
. “ D ОШ Р АГ р nw unm. eir Re MBE int 4M Corus АТ BULL ILL | ез, RM ra Be D је ИСТУ Мера 
1 < 4 7 ГИ > eh ped oa 1 А ТИЕР ТЫ БАГА „э фа e” Алу gala: та (резени + la АИ РЦ РЦ РА t 
u ` k ; П ТУЫ о “.. ТАЛ ЛЕСКА КА РУТ У УК A мари diens бүрү A E Е Polen Fe OE 
E " У E ) v n PPS eo oe А А За айч. $e. MAS ule. ЕЕ ЕДИ АЈА а СЗО етте ЖЕ PX" md ее T 
' га т р . О РА TO k TC SIT Ym a ЧЕ ТК imp гази ги УЕ: Сг Pe an es Eee .. 
б, T ИГ ao. a ñ NX TCU NT | 2 Te ee РАНИ АЕ ОД ЗИДА НА У РИЧИ С dnd xu 659 eN pite Ü bua. 
О + ] tt ТУРУТ в ПРИ ОТ ete i Ee ripe mrs Prey (P HTC Белаяр дърт Ы ааа ын P аа атар 
. \ " ә А А i AED Wa m: w a. we. LANE Г ^ АТА С a лді. 9 аа РГ е еа ~ очурсо 
T lo ter A A dee a. B рту Mo de ИЗ ЖЕРГ ең d aree 
О a BaP Bl Dee FOL АСЕ ОРС EN ЗГТ Ba Wp She А ple ete МР p MM кырый т Дн и zi. s 
а .. г H E ch гү eer t ELE ee brary Ди LER: uU Lr я о ад Apis UM м y 
"LED T а Је РСЕ ЕТТТ t> AAA AAA А AAA A MA Nomcn. ve P Rmi cr А nid 9) Je a Men ee f. end 
та . : REE ua a hem A E T T e AM CPP Y VES Зы ga ча pote uo dtr · лы 
AME. . ліз а 3 ома часа 31 ма фа д ала Аре, PvE [OLI РА qn EI 4м, rh ааа 
9 safe ы " AA m POPP Ад УМА, ETT M P ees TA Feder, E: Бон by tp eee р 9. Манао ы Бае 
ОЕ D t П П 1 а ЧУ КАГУ T, TR id PM rv TI 1d Thr cam гы e Қы өз Me 1а piri LA. edema ТР И чи 
ao n B " - " 47 МА», ак» АУЛА S OO АГСАН it^: тарамі. О ida 
Ы T^ "An г А АГ а ПР РАНА рече руди S E ITA ARA A ы а И d bien У lb 
.. П EE 1 . е n "а AAA РИТУ 1 Jot rire ale rt EIN Sa rata Ч 
р Не до . Е '. М ru 2 y AA CAT ee? а ФИ ум КА LP ASS РИКО ч (ira 
.. "E D А я Р; ЖУРЕ РТИ РТРС РРР Afer Oi Ce LC e he LU EE СЕ жаға, ся 
. a О ` ва x €“ w e. E t t4 REC ME E КТ па БТ РОКУ A к. КАТЯ їл дь 6. MATE СА Па Hs LA 
О . at ПЫ ПСЕ D 0 AMI е n «~ њем Y ~ ARAS a We et hs РГБ ОАА e 0) жаса, bh Дари теле dj Аны тте нын; 2n Рагу 
О О E ү га! SN Y Рим дз LIE, EN DOCET EST о У ОЉА У ЛАО АС aE Ми PE ааа, Ва! 75 ЕГ И ед строени Зы Únete 
1501 Ј О ` LT ILLUD р PES ID а & Bache x! О РТ m Iuda LEER ү ү А 5 oU LI dais 
[] В П П Ц n 4 А ELEA 4 D] 4 TE Tn "ім pis $ А ia г EI Аласа айы LU от лама E А үсіевм ж” ze PIA ка vea y ink iei а а es 
.а О А Ш t , ПМУ CTN A "Ы ACE Ca hy Ory Peer ey are ewe АВАРИИ ^ KORA iq умит с тале, аб долот обо вайо тә ума A ` 
р ' О А b ‘ CN Ы NES OE EE WC Т АУЫРЫ P Да eae west ee rrr one ЕЕ H SO. HA к сен ААА А МА атаа = арачы қы н АП meda er 
y O Й O i те. E 4 ; DU TV АУДА ы ТИ re К i НА El + и 
Ы D 1 i "1 » В mm "mue ТИЧУ reli Wt СИРА А У тфа Ат е yir m уа neh ува А „Ам 
d « De П М ^ E [TENET 149,1 154 tote he ` П ЛЕГЕР РЕ РУУГ мо ај, БУ] а ede Sain À pe ЕУ рт Сл и паре $. 
Й B + *« 1 Г o П в > 22% Vike ai DT" "m DN УЗВ Ys er a тте > a. jh, arri е ж ДЕ м. E 
u ^u ' ғ... Ы TEE uf ‘ase K s б ЊЕ ЗОНИ УТ РА У ТИНА ЕВА АТА ЛЫҒЫ yv gu en Loi ri rta db «> АДИ L би = ri 
О G LES * 4 "he P О Ы е ао 5А й h oa Di Pk ARO а тфа. ЕЛ грее ра А EA ~ к ај БАК „бр а Таљ m OPEP: 
A Й a PATI Ta E мора фр Ра Ке, Мо ВИ 6, мл ж че ИПК Chere er oe үү күлө у СИ еи EE V wer n к. къ 5 25 
= Џ . 441 E ПЫ р ELE VE. па а ма 4 DET TYPO e ДО н ТЕТ ПРУ wes toS vna БЕТАНИ Т T ues A КО ode Т Ti Ha e 
D E .. H Ts d á | л наў гм "TR я: леа НЕ P" [os As Abana) ra) [NP ч слуга мо А ммм mi ANTI раме чира ai йар. У PU = a 
П M L mo .. ЫЙ, LIE E E - ұйы ы "қы 7856-9406” Y e ads Aj Mo: УСА упор НЕ EPI ае bayer > vy LL = иа 
DET , Пт» UNE BET М ag Tae bas * "LL DECRE y . Дъ “ CAP NON AE LS Mt eS 144 T f a s лэ 1а ды e ss Di SL Aito dod e 
D Й ' га ы; ; • mm ^ ELE PO 4 47%, METTE A RE ED d дари млада D ту ТАА rng. ice ыма жалы 
Сү yoii ü E A DLE СИ ДА SU ets eya, у руе МАЗ ҚАРТЫ озера: eas Е, DT “мз T Е9 me Las eed uei 
B E " о . . ee TA IRAN TES СЫРЫ 3 Le катрае Део уља ab «> з arme Re^! ns xdi rur 
А D р E m ED А as . РГ га ледът а г» РЕС BO ы ЗА РУМА Che EE PEE» Ада аа pet oe cy petis Panic ЕТ Ra 
L] z . . О LE А, * » od NTT WLT РТ КА a " а ч" a eae at Цели ETT ВЕЕТ at Hows iM Po a ААА » Армат ақа rh mm 
М s uy T oui р 44 A EL ЛИЛЕ ЛА ЫА? зу" nt he ч Ye VERTS ines LIU Pr TX d "n ГЕО жазала pra фы м! R raps 
Ш О 4 4 M. LI J D *> v ed арды. A pter PAS. y AUS I SAT a um IIT T TX eee PETET ears Е ии сы aida гм e 
' a ` 5 у ' О А аа таға мон так лыс OR кесе к s дим оу мм им ао РЕНЕ УЕ РО ы ых 
E ' А Т Р ы Ш М :.:. У ~ 2 236.05 АИ и ХХ P TES Да МУ o 5 дан he usará М. 
. " ja sh pri’ Н ^ РРА О р . ДА е РА А | 5% reer ean) е "UT КА ENS X ама?" ри tA. чу y ат да ЕА 
ЕЭ ' 1 И ы я у ы. OM id Фа ВЕ Ат rat pa más кісіде DIM veas SEES еее rh 
. а ы .. E t 5 sa í ~ a E 
ОТЫ р Ж: ы” £X he AA ка мал чи . AA db) ET m ds, a EAS че Zodiac ha da GA ihre 7; D За 
"оо, ^ a A | 4" s .. 77 « '6' „4 año, ds Pede УУ ма ба. ЗНА ber E УМ PC LP rimane IET W idi Ва ЕА и ады 3 пра n eh inb keen 
. .a Ys ще с ПОЋИ ~ Е PC LP ' “715 ACA id % Бе ЗУ АСА чо ime, 1м Е fr nhe deir uber с трдн га 
Й 1 1 a age 5 ПЕТУ g Fi PE NP LE T pA M, ANE ии y HE vro ра e “.х.% d Cl a ak ЖА губи ағарар 
ШЫ " +. y " аза Жаа: e 0 Й = *,520/X 425 Ww 1 ЕК СТАТИ Г ES Y S n Ж. e. + br ДА Хе Ж үс, b p 
g А б E TE x; тм р, FL. Aa Yo rnd MOTA САЛАНЫ о ем У ai Ны Аы, X Са р: 15 TAN LE а 
TL "er Ure n И PL ки жр ` q ер 1 Г? (EM. ТЕ те um. Po a ARMATA wie ue o ya t MI motn ES et Де M LUTTE Ат) Td үт 
Е и J ы E X УЗ У РРА ЖАН ned V iu АИ d e: ei d AAA A E "hes RU 
ЖЫТЫ ЕКЕ кунү ЖУА ала ее ni CHE CR а н Ki р ДВ E DM < 
И С PLA уд E а ии ба АДЫ Qu Бе Ux 
ТЕРА ЛАНЫ МЫ. "oor opu P а гара Ы магог A Af e ATAS AA = БАГУ, aa 
а a ымы талға mM EN MY ТАРА Xd Т АСРЫ ОЧ pop nu Г Сы ы Кы Fiber Ур m фы NEL 
ПУР 9 TTE СИНА ELENA ert А а НАЦ СЕ r ае EIL RE he д а m кыа а ^ 
O CS ES i cry AS щата печат Е ШЕКА т oboe: АҢ, ОУ Jis hj 
Py ат NET: wu wae te FII Ma hye P y x EE на 
АЕ NET NE PL Sue = БА лаа ИЕ Ha are я (+ ЕРУ www 4 i mpg S S Mire nui "e= 
и р к е» LT De x e xi E y D АУЫ D ПАДА УКЫР УЛДЫН E RT. 
E Ы di Ez] mr a t an 
AD se EPOD a Uy Pr e pre rst ще қамы s AREE 
pe DAN bs T ~ ЧУН E E EH КЕДИР CERIS ПРЕКИ Ан de 
4% Ру АНА ‘yr ADA “А саны "ou" yA MEC иде bes dedil 
"I ху LM RENE Ж 333 РЕДИ ИЯ ЧАЯ rn EPOR TR) қаса be rey 
~ О ths Y » ti 24 ЕУ ar ond E ERE д № acpi even ~ Lt her езер do 
ы В 


А E EA AS 
AA DC Ede CEDE Ga ta РН 


A tad ie 2% "bb n Це 

















































































































































































































































































































































































e s ча a DP 
р Y - е а лы ааа а e ARENS ЕУ мм чч НАИ турс reni ро 
4... FORTC eel NA ы b or М re a h NE E H: О E oie d ЕЕ 
s E 1 ARI TAS е4 МИР РА а. Ж: 
E zi в H p г bi A al aya ha dl ^S HIM күз Bo tO EE ои КЕ Г се фид с 
- 4 Aaya? ú w tat li са in IT ER T M CR RE abd br Tin I PLE те ы ДА аучы 
ГА мо 7 Рабл МА АЕ и 19-4 "EN ЊЕ MN ti} ЕЯ 74 ча ЧЕ Le м rhe 
г А НУУ аруа; б n NES аси т: £u а Hert: 
4 АЧ қ-а Фа Ak tE б H Ало OT М үнү 5212 errem A 55 4 PEA >i ~ 
u ко M # m ДИРИ d “Хғ: E ^^ КУ 4 y эмы 12 У b Ett VELIS 
] g СИ NP DD EIS Larra А p А 
Й Па # ui A 4 '› i A Lone ^. 44 fens E 
А ПОСТ E y таг MAL 4c ed, П y ~ "> АЧА У - 
= ы sf ft "^e r ed "d i 
СІ ГІТ А d LY dcs 
44 %ы, y 
m. 57; D SE TW Е да Е RD Hn 
Ы * < A PAE O =; 
П n za x" ан РУ 
9.- 23 4 ^ ъч р У НЕЕ RR DER 
. Ру . б а. M ж e HT с T 
MUNDI УЗ [ES puta К TX pt i ф ҮҮ НА e e A PP CE 2 Rt PAN 
ED à ОВИСИ ` E . P . T" ПО бесно т На У а че A dS пай n> EATEN Ro FP, Ее 
r . Й Й А Ы Dan ' k r; . E у 27, rq "K PET “22 atas Pu, (ET TR 
. 1 А ` А s 4 47% "yon as P .. , PCS Ре Р * хуу R TIVI ДЫ. me" X Ff Н P" 
са s, s © a DC P А Е И ТР Ба) y tal tag fet Fates ig gras EY НИ Е pao 5 сыш 
ED 0 a ва , " E , me 4 А 4 А ПЕ . 5 DAI mr tu po DT Р 2 Ы бе a ao, Е еы ыч 
E PL ха ТЫ) ь . M. П + Ln "n E қ О РОА RJ 4% EL ПИСАО и м ..; 7% 46, sten? ки 5 y EU ЕЛ У E TA 22 £ ра 
Пе n qa ee p n 3 4 ң tup ОЕ "m РАНИ: ive E Ир И PA А: FETA At га Неда бір Pi FER 
I E С I uh MES REO Sr ER 
Eos d у 1 k < E^ Zu p AS Fur МЕД A Fo AE un rta AR T e 2 ж. 
СИ ва ... an HL D depu s НР $$ rl Уф: p bJ x» 
a n s. s ЈИ Да А Tl de ЖАУ. Pe ee AS Ee, P 
. m О E ИТ. у > ГҮ, Imre ET © 4. Ди Ж ЕР Ey. ИН # 
) А ы t / + 5 T: 4 : А -4 
ar. Д с Ба AE T ӘР 
О ü , у PITE Y iMt TE Dor rH А ра ТЫ ге Ме“ уа таме ие 
. А . МАН тар ЕЈ қ UP ер ME 
rho ^w . А MET "T *y um IO ES Vo E 44094 2 ese"? 
axe .. б ма СА + Фа n OE тека P Ç E Еј ар б AA 
У mou | Жек; НЕЕ ак: Pe Рат 9. MARRE 
I F “ F ү ма ве а Сіз LET Fi ж Vu oem EI 
i ји. | ТАЖ SAT LEER B C TAS AT DES ERE EM fs j VES: ТБ. REDI lee 
n А ан 1 INA TO Rtg tt теч UNT LENT AE ACE SEL s me "ETE 4 opi я RARA A+ ESTA 94 ыр е d 
ШЕЕ Д ` у d г ғар Orr РУФ t, а ое Г er “ы n к даи Мона! ES RI МАЊУ tr pens | 
А P " ' ү Du е rig ee as Ми X + f amm ee ГЕ Y РТ ts He HORN EDU 5 PI PEA ci RE Mir 3 
10 ГІ NIE 4,4 Й Ты = Й 1,5 LL и: =t THEM + т" АРОН Iu, T " Де AUS ~ г НР i 
. doc .. ` y Ben oe >» ГА E ДЫТ, PN Y 4 `, за E с C Емен rena ХЕ, "Eli ж. AS Res а: ен Е СЕР 
П О О АЕ "ee [TI % w Ir. APO Tr a Р. ПИЕ + j'y РТА P jx d Les de nr LE 
А 1,4 Wey 5 В D (7 yore ee Ë 
Й . "E t ' PE Jar PE C Er аў bd ey Е из 2 РЕДА АА " ОТТУ Pd d bed e Tp m Pu au 
I ' r4 are eS O А Ex “ag Й 1 Пи 9] М? „0, 2" 42 ач кӛзі [a ші; БР МЕ К v an qui » ЖҰ а LN Барје | 
Й .. a4 Aer Й s TO ТТЕ Й eo PEO O at et Ne me m oe " та 4. Ы; я ЩЕ: ди. 
Н D ел £0 Й ГИ] т — 4 ae " "ADEM DID РВИ. { igen LQZ t MATE EE я vie 9 ro." 
no. D š. Г а оо 22” ГЕНШЕ ТГ” 7; . 4.44. "TUE ПТИ: e sr, e. aoe и. w Vh at turtur 
sa Р g баз . 1 Г " Й БКИ а ‚б, ET A ИП Leese ma ts СИ EN m ыг] т dium 7 ки; 
E ^ р жаа 3 UC CE $2.5 91 г ДВЕ И РИТУ ва Deu ТАКТАЛ ШЕЛ, ОИ а р L 
р Ц "D : mn . n ГИ) . n Weary a) " DNE Tr as РЕ 4, TS "ru " р E "an rà КИИС ШШ, TECTO "ECCE LT ACTA НАД fit". 
А » H КАП " ОО ов E CULMEN" Eun П ud я "PET ОХ Г ер 4 б .. ИЕ MEC aco sts у 
g Ы n 4 ЖОЛГА ТШС ДҮ) " П .2 МИР e ТИЕ И ЛИРИ mn Qe dee СК ТЕТ КУШ E He PUE A LRL А $ БЕТЕ» ` тб ерген"! е - 
.. . 4 П “э үг Во РО то A "un ..! И . .t qa А ^ .2 t EU АГ LE TII 444,6. КІ 2 aaa Кш В Wed po фа A PT P RAM f. ay тър: 4 
ae б ИИИ О EN а И ... DEC Pe rer pace Е ЫАЛ ҮЛ, м M LT, G ETEN K DP M SALTERAS Ey prs А а, — РЕМ 5 
. . “и ... "an t >» PKT AOS OS O пог те aw y + Тар ғ сан . СОРИ ЧЕ БРИ] Ft en FE pe +95: 4 Ра i учи 2r. у iE LES 
СИР га $1.87 à et. "un ба nmn . [T] и 6 re ОТТИ! ht ПА в T ҮЙ x 4 Ld M. ut dm T AC а a ber ac Ы Pr EN SE 
н ' + П ПА П "аа П "m 4 А Г ПАРИ . Lu *» ме NR у MAN Po TUI Du id LET б "i Дега Р o шры ot hed 
Й ..a асо ИО Brine (52 E f Р оа a L n AA Тее . “2 P и и t EH 25” ҚАЛАН era ona brio 
. Й а. ' O oon? Pe D .. ШГЕК. LI ТОС .. 4 E I gue BAe ol eee at zr ДИ pw rx i Pi EIE dB Bat d бер = 5 Eun 
Bu П boa ги "n . P oe E ra yo "а "Jab od a of aS 37; ж a y... ee ПАК) TUR МЕ (ieee, š ге у Диви aar A Үдер 
.. ГІН! . . 5539 LI DTE AT OT EE И DEALS и МТТ п ГАРИ TII Juv ro Pim tet ГЫ АП ЫЛ . E rd en И volets rap Pr dodo d E Gi 
LT ПЕС) Па а а Cr DE Ирис ва » Ти ТИКТЕ И: 4 Uae Ee у ° ШУУ 2 СФ ИДИ У D 9 ВРЕДА "mE = PALS BRS. A LO ¿E Ea aie АЕЦ ыы. #8. D 
ЖЕГЕ ЕТО OA E "ea, 2” " Ре ва И ОСА sus TA O ШЕ олы s Pew ICO quae РАИ нету ГЕЙ 2; "m Зем ТЫ Fo жу ыч: RU AME I pL e E RP eni Per eme 
в. P КИ. ED "OS P ІЛ ТЕГ NINE" Г iut nd ПОДА ПЛИК ЕЕ КИ TASTE ТР КДА А pm 0 а > pe ET. 
EN 5 > ҚЫС. ^ | A м го зол imo oc "n A Tis Fe БУ ле I БЕНИ НЕ s wo EE E PU" PEL ay) Ка NA Е эн yea y ES ol aate . pr 
- a АКЕ, Hn NP Й Pure TT «а РИ в уши edo 7 П По Bae HF te sa gohi aY poh Ма ја nn er ні EST A ” TZ P nated ра Fa lay edad al ct dc Да hi ui ul D nee 
М L| wa O IO КЛ ВИ В ey [EP а Mee dae wt О . да “а DO A O ОНЕ IA IR O ТЕТІ ТЫ TA маи а, "ч IDEM FLAS of pt eee niet saat TZ a ia тл; NEM A 
П СИР ДРИНА П P О = СИМ О ' , АДИ Й веки CE) СИМ t ђе Ae e рай дави La tye Ft rh “іі... 1-0 IL И FAY P» rore t Imi > FPES rl Cu 
.. А С Е ама L " oe P mn (th, MP СМУ: раја . |. РЫЛАТЫН ТЫЛАР UTE e be he ee ae nar ed асч ана ра 
' ғ Pe a +, Ú ED ; ` ОНО ИН ИГРУ АРИС ОР ШАР en LC sii TP EN] Ped. dd ii ИЖ, ите yo Mes АИ ш 
П P PETI П О ... ñA я n DOT ES то су тори отн и Па ІТ Yn aq» wt pie an РО а И а o ak 
“7% TE. n TP 17. P 5 ny на) rss 219 be "m г + Нух и “i be ree War A ее 
. ПО PE i s па Ч en Tay ЈЕ у Р АГ wy ygt ~ a” = 
| ое | 2 6 Е 
ее ç АТЫСЫ; ҰРАР MSN, МАК А кр Pp ГУ» D H E EA ИИ Зн р БЕЛА p" 
Led мат а. > р И жат Le КА Ки re Y E Fe: he By ро Quo ti it 
.. .. 01 РРА War ON | ЗПР ДУ РР 13и ТРИ pr a Ау Ана о А ry ver vend 
MU 5 P n DP ПОТИСКА ТАЛЕ КЕ РРА tr x pw T. Н ещ wk пи pp à me 
П LED За 1. ПЕ ПАЛИ ТРАКИ NER ТРИ Нера n "7 Us яу 27 О M ^ раз ig go 
P " d m ПОЛОВ ИА РИ Шы y TEO VIP тати E ИИ а y £ КЕ уљ поје PE C gin 
E ГЕ: Й ET ani y 41 vi PA 9* od P^ И S A EDUC g" Fa y D A i ay Ша 77 ыы Заро аа ALA е Par арени: В) Los 
| Й . ЕТО ТСТ MAS ñ 4 paa 4,494 400,7 а; 9000 03,0) ius M UTE E AP БЫ n И PSU КИ КАС ти regu в 
б "uw Lo E ПУЛУ РА У УСКИ НА nie во Дони: yt "or е nra P ТІГЕДІ) КР НЕХ аи Leg teo m ec 
HE 8. . . O CS LM oe De ИМА Bees Uu ОДА АГИОС АЛЛУ НАЛОГЕ: DT ИНЫЕ v И 1ч р, а пи: р Ту 
EL i à О СИ) NOCTES DEFINIDO IU ILE М ЖАНР АЛА, dic de ras ДА На © regt n Б УДАРИ abat 54% ыы ја аи {д wrap ДЫ 
А , . ў Pa ТЕУ ТАТК ЛИГИ КА: VEO COEM ТРЕ Г Н ы. МЕ poc E E n cree qute ud A ee FA Ei m, арх Emu 
zi 4 LL Џ M , КИИ СЕТИ DEI LEES NELLE D COR O F P yo aoe Ж О МАРИ ЛИГИ - ПТ) DIT 5 DW Tn ETE LAT] PI i DAR. аа eos р M СЛ я Ce sai od ai 
n t n ' П * э„а в E | А АЕ NEP" ана ара Тва бүз 4 цев |, НА ИР; ЩЕ ЛЕО ри АДЕ, ae Р ке Td ЖУ Амен ннн vx. Др. Nah paced "ic 2 PI Es ыы m оқыс ча Ll e де 
ate "nm . . ой 1 moo "Om ПИТ E за АА одаи е Vie Spee 4 ОКА Мг СЕКЕ ; E i ен У EN CAY A i ld Ink сы 
rea ier 4% 7 Е л. a) Chay pee ES SO РЕ НЕНТ CLE теби D aM TO КОНЕН SA E at Pic pre Феб; и При Ack ера 
. А Юч П sfai n. tio п , "PP: UA" ANT ПРУТ Р Пица аны E г" Аы өсін iade aao 
П П П | » T | b, ДЕСИ Ú D a 18. eo I PE ELE E d Mee E ET 5 4730224 INL рид alcoi ан 
' р Б К. Фо E E СИУ RECO УМ ПЕТА а CERET] ennt r peu БЫ 4% а. ога а у: en есеге 
Й e " " П EI в..)0шв.40%.) 4. Й M D m ЖИ: |o РАТТЫ Лаже "Pm : DU АУТА e unie и е 
Lxx P LEES а A NE AS A EHE Aer AIDE и 
0 М L р 0 в . а О í . ж“ ТА; m H ath H X пёв?у д А Ай hay ve 
А А . 4 ДЕЕ" РРА Не la - nom mbr) 1: НХ 4 VIA К TEE АЗ 6 P ec Meer a y о ДА АМ СИН Ми Масон әр Зе 
E А А ui ew os 2 y >) a a o gogo ууу ЕР RD re ee th yy Т; 7% od = pip Mi ir arid beri A 
.. 4 A T (4 | 1.4 + + ET sph qe oe 3 ми "J Ше T Cie желчь pr gos ИТД тан ЕРА СЕ ТАЈ е. 
.. П ' Й E И га ПИЯ Е ows ot raphe pe ari per * Кун ГАРА ais Ear] des DUST pm Қасар Tr Dm 1 7-р 1 
„э О уа 1 y + Г 2m PP n E РА ЛЕД, ут. pt pte per vie ago C ДИР И иа: "ns гука le D M dl o tmi Pat 
$ "n "n С n £ at DEO г. As AD И А ИН АИ А ИЯ Pu 21.20. T ipe pic С ње 
rs y o ND г " | EMO Mer Te T. PL Cee AVP eat’ а Ue y cu yt at У ест rea I" мнение DIM Ж РИТА Я АА oa әр а 
4 PED б яға баны 59 И T ME AA E ES SS СА АНА АЛИ РА НИЯ . РЕЉА ВА УСКО an e i x pes 
А 4 Ж. СЕ“ EG Pa | ve qf MAI NY ИТА ИЕ E P Nd np IES ЕЕ ret LE m той аа IV AP UH pote d A Mr E yup UN e € mies hes Ai yum 
МА] P Й a ES IN "a. — a aX ЖИА, Т eee B® OF KO See ми СИЕ ТЕЗИ ШАЉУ "^ ey CONTRA EL db SOR DA или г" ПИ дек А sa 244 pori eto 
а а ПО оми .м. „ i: . "v ИЕ. GRECO E AAA R lod и МА И MEN RT ES E AL n rm иа АЙ nt EAE H po cai: Lans 
П ‚у , t, 4. ге 92 1% фи %% К O AS ҮШ RC. а JC IO DO ИП ОК АТА LA P, АЯ PT la я Solow 
. 4 b 0% unn ШАР. D П РС 9.4.4.4) әфігәі? a f£ a8 og di 7TÀ O По E М И ВОТ Л VELA; mir LI ИГ) б ПА КИ Ани 
А í. " f И nr en % 9 " SANI ya t РЕТКИ TEN" RIT И жү PNE Dario. raid ds m P T ua ТАИ Rm E Pier r и giis pate лы Жы! ç 
А Й А " i H te Га е. aes 84 E i ei ir ek ere he | МИР MI 4741 ШАМЫ луына ры utm iul har fopra da! [2d А ҚАТА e a 
E E " E Ке ç 9 x Te И. PE ә” аб, ( TUO д. м da 1. E .» 14% Los HOT IPTE LS nn Pur ТЫ ам os DT ее Pes Tet as ih asp pakka СТ 
TD ... y EDO TO ON Gc A MAE LULA BUDE ead АИИ АИИ ИА ДАТИ а И У НЕА ЦД Mes БАЧА nm ГА ЕИ Андор эй ae e 
n а и Om See e ç` m А OL. ку, ЖОЕТ ИИГЕ РНЕ LTU Lea. Pont? ae НАДН ОДА ПА А DIT TIES CTE ma peccet iii 
Й Orem y Г ПЕ А в г, s "s ET ТЕ LEE ORT] Й OE y LONE " ES Л КЫА st TM y 1, bah А Ы АЛЫ S ЕДЕ ГИ "v База PELLE: my par. hj. 
а OS ы 1 on” з! qe 09 E ETA O O A II ЕТ Ag fole ” аи И A ишел gt i ағын тым Бас азу ді pepe M aber 
" vota Dosis + TEC "uS бг, WP YE ие ми Ц ПИ G S Kp y, V y'en I al eaa ri Z eL d Lone aM TE PE RM i ue AL pale ULT eie bos БЫ ru Aaa A 
. ‚ә Са и ИГҮ А ae eo SS TS ЭСИ TEE па ЕДИ A ss БҰҒА he. mon Leia o quia eo. E ез pe eS 5 
D E г» А 1004) $ “эл, МЕ XE рос би №, №4 Зеро: КЛ мие И prb d Mir ру m ИТЕУ АЗАҚ үру wey аи АЙ, но h temp АЫ та ne 
rr org t Й А 1% ' Bo. 04 пу Гинко пе они AAA A Gola how ТЕКТІ ЖАТ), У уу vw сила Ше МЕ) прати ока ¿A а) ње] о. а 3 
О ' n LI П E 0 е} П VUE ЛУУК. ee AO le д ОДЗ DERE тама: ДЕМИ i изу ыы a d дуо. КУ M ве. TI КДА ТАТ АЛАМАН МА а Да 
1 Б Й ' I Sa А $8911 1 9 95 4 Pr "и “УЖ ТАЛЛ: РОМСКИ КУ С ЗОРА bes на m wi ДИ Пал Er. ОРТ My ee ipli cota pb sia ii UN 
n Й А 3: 3:9 ү ТІНЕ" "па. > ГО à ОЛА, w ` ue 34. "n e Lb cde ME "AJ an чу РОДИ ате Че TT рН Lo ad = 
П " У Й $e PD П D р T NEA ИЕЛИК ТАГ АНД» TA AA A ad IO By o a pe rd d => 
PEL x c я , ТО. #11 По n Ай ар я ек ОКК Инн “%; ae и В. ка ма ту 4 МАУ Роа Ја РА КОМА е " поље правди ин 
u " i . rae D | . m ʻa E t i "oy 4* t YN AL L ДИ ж/е; НОУ РИСТА БАЛАР A ОНА O Ыы, g pu TE Po mI E ка ава 
D УЫ” m [NN NC "s ЖИА а "i DET) В СЪР 7 АЈ pr Ги Па ts d МАСА АТАРУ А ДО ДОСТА ТАКСИ НОВ БЕЗ ЫЖ А41 DA сдан боган 
D : " ° 1 «а... фа Жа, i LEBEN АЛЛ БЕГА] one Pu НАСА а АНЯ т 
LT . И к о m LECT] л E CT НИСИ А САНУ АГ А). DD Q Mm HD YE U] pis: БЕЛДИ a br GOO EIA АС OR E E qlo mse, dei кезеген 
n 9 E E Й Й Fa " n i " d ПШ y ПАЗИ Та МИТ ТЕ) HE DLE ИДР. ЖАШЫ | WKN Cory và gy Od не 2 АДН. ад. і 
А $a! Е? г М ыға? уде P | eh ТЫС аА: е CTF ЖК EA FLED CLT CS CURE Nyt. ТЫНЫС ыны әд а, v 
4 P 2 39. 9. 2, _ ГАРА ИСТО ИТ оу TR | (GA реда LINT Ме Зее LE АБТ! ro Ardo CL rl E e 
С А š 4 4 Pir ТО ККК КОЖА ИНН АҚЫ БЖ АТАСЫ Дукла Mol HE. эген n А-да ремете vad 
. . D Й 1 П sae gan ek UU] еи AIN ЕН РА ЕКА: ПЕТА Mei ede | pre eM мы 
. 1 T + А a w ATTE ISE OA A A A UA AA a a К Еу Sotto СА аА hbi арч eR yak: Lorie 
.. ' ' г "Er E | Ў г UE тр” ROTE d «он ен АНИ ДИ АЈА ч Трн гази AAA ро р рат 
' ç a e PY м Ж ET A ДА ПЕЕ И а г al gen goat ue 
' i 3 у 5 HS ` el Ki E ri P ; ANTT MP УДА aE eS oe whey реди E pu Ы aay пет иа ИА pue 
ULM. MOLESTI MEM T Ажар га туп А A БР Аб TPUT š р, f MAT Wir ri dme: НИНЫ ГА TET HE МИМИ ИА e да А > чу Заб D opm 
| DNE ME. LEE DURUM pr CS ЩЕ ADU н уе оаа АСЕ 
и ГЕ: 0С) EN ul Ln y. i л ` i чат, Zn 1^" АЫ: 4-5 ie 
A a э ET M - Ç ria rn P y е du qr е Фро ДА, Аы; DAS ‘out ov TT ы b ГАРА а аа DULCE ТИЛҮҮ ПАРИ А расата 
: М = Џ ү I i ГА TU C ПРИВОДИ ТИМ И БОРЕ ИО iret Tm) RY ie eem m yr А НА LU na ПР 
КИШ; Гия [En па а + DD MO y Ол [E ТАЗА ТА Га: PIT RC aT e U Дже АЕ, Д жеме 
р { ^ t оо м D MINA AMA Dd del A г ПИК а фа а +: TRT ме Ки DS ke th Le ie bio) c Е ^ 
А ` í Ld JE CPP E $ i ol TELIA AAN, Да | EA vá ra па аи Да ея f 
SEA ^» "f, ^ PRORA ra FY Teg ка P КА ИС КЕЛДК УНИ АЯ и TAGES Қатар кашу 
7 ай d Ut А.Х е а и Lah A den SAA НЕА ШУ даа иу Пака e aa ылга 
р ч Р. П . . бе, » n3 мд ya 7% Й пса KA Г ым NU nd E a eee то 
р 2 и 4 ҮҮ EN Тай, г НЫ аня H ҚАРАЛЫ Ay kakak) дар але Ка иу, pate ЫА пещи T 
nw А ? ACA os ARAY Ina "91 ty ae ar LH o d ZR BY hea x ARS eT "E eT 
vo А Ie a, Даљи а МАО e КЕ >ы MUS AOV и а 
Й W Й ., „пору ПИ АИ у КА Цао ОА мила en 
А П ñ n $1 m Е: D du" жә LAE, I MA f, T? niu rn Он ОВЕ: + T "i a 
у ы О ; s EIN а! АЈА, we LES | hc A g тая БАЈЕ и А ce 
О ` " б > А P Jl C q^ oa У) "poeni Hs ЫЫ v v: y roja КыЗ p 
[EL D П 3 р қы ә! IPE А Sa ТЕ UE 94 TA Eur DRA 
р l AN NIMES P de. 7 8 V, 4) ж м” EH deg VÀ LL 4 по Te" eon ина МАХАЛИ 
9 уй А, Ld eo Ww WOES Vor n FE £x A ee теқ ЖЫР PRA қт Ы АМ ed ЧТ ЧТ А 





a bY BNOA LIBWAWY 
V^, POST. irs UJ ATE SCHOOL 
NIBLIY, veLIFOnNNTA 95943-5002 

















NAVAL POSTGRADUATE SCHOOL 
Monterey , California 








THESIS 


A STUDY OF SECOND AND THIRD ORDER 
MODELS FOR THE TRACKING 
SUBSYSTEM OF A RADAR GUIDED MISSILE 


by 


John W. Williams 


June 1988 





Thesis Advisor: 


Approved for public release; distribution is unlimited 


1229220 





| UNCLASSIFIED 
ICURITY CLASSIFICATION OF THIS PAGE 


к 
Í 


| REPORT DOCUMENTATION PAGE 


| 


IS REPORT SECURITY CLASSIFICATION 1b RESTRICTIVE MARKINGS 
JNCLASSIFIED 


5 SECURITY CLASSIFICATION AUTHORITY 3 DISTRIBUTION AVAILABILITY OF REPORT 
| Approved for public release; 
b DECLASSIFICATION ; DOWNGRADING SCHEDULE ЕН 2-21 


| 


| PERFORMING ORGANIZATION REPORT NUMBER(S) 9 MONITORING ORGANIZATION REPORT NUMBER(S) 


le. NAME OF PERFORMING ORGANIZATION 65" OFr CE SYMBOL 7a NAME OF MONITORING ORGANIZATION 


(If applicable) 
62 






Naval Postgraduate School 





E Postgraduate School 









75. ADDRESS (City, State, and ZIP Code) 
Monterey, California 93943-5000 





ic ADDRESS (City, State, and ZIP Code) 
Monterey, California 93943-5000 






a. NAME OF FUNDING / SPONSORING 8b OFFICE SYMBOL 9. PROCUREMENT INSTRUMENT IOENTIFICATION NUMBER 
ORGANIZATION (If applicable) 








E ADDRESS (City, State, and ZIP Code) 10 SOURCE OF FUNDING NUMBERS 


PROGRAM PROJECT TASK WORK UNIT 
ELEMENT NO NO. NO ACCESSION NO. 


1. TITLE (Include Security Classification) А STUDY OF SECOND AND THIRD ORDER MODELS FOR THE 
TRACKING SUBSYSTEM OF A RADAR GUIDED MISSILE 








2 PERSONAL AUTHOR(S) 
NILLIAMS, John W. 


3a. TYPE OF REPORT 13b. TIME COVERED 14 DATE OF REPORT (Year, Month, Day) |15 PAGE COUNT 
Master's Thesis FROM TO 1988 June 139 


16. SUPPLEMENTARY NOTATIONT he views expressed in this thesis are those of the author 
and do not reflect the official policy or position of the Department of 
defense of the U.S. Government. 

7 BOSATI CODES 18 SUBJECT TERMS (Continue on reverse if necessary and identify by block number) 


FIELD SUB-GROUP radar guided missile, tracking, missile modeling, 
—— | ч COrliOllis equation, Kalman filter, aircraft 


modeling 
9 ABSTRACT (Continue on reverse if necessary and identify by block number) 
Bis thesis is a study of missile and target parameters used in second and 


third order modeling of the tracking subsystem used in radar guided missiles. 
suidance methods are analyzed to determine which method is Optimum in a 
search for an "ideal" missile. Target parameters which have an effect on the | 
nissile tracking system are analyzed and a target acceleration probability 
nodel is discussed. A two dimensional third order tracking model is simu- 
mated utilizing a Kalman filter for target parameter estimation and 
orediction. Linear second and third order tracking models are simulated and 
compared with the third order Kalman filter tracker. This thesis concludes 
that a proportional navigation guidance method, with a non linear third order 
tracking Kalman filter, is the better model. Benefits of using a non linear 
third order Kalman filter may not override the cost and complexity of 
implementation of the model. 


20 DISTRIBUTION / AVAILABILITY OF ABSTRACT 21. ABSTRACT SECURITY CLASSIFICATION | 
X UNCLASSIFIED/UNLIMITED O sa ES APT C DTIC USERS UNCLASSIFIED 
22a(NAMIE QX RESPOKSIBLE INEKUIDDUER I ` ` 22b. TELEPHONE (Include Area Code) | :2c OFFICE SYMBOL | 
Lc LED 408-646-2560 _62Ts 


Нагота 57 










20 FORM 1473, 84 MAR 83 APR edition may be used until exhausted. SECURITY CLASSIFICATION OF THIS PAGE — 
All other editions are obsolete ft U.S. Government Printing Office: 1986—606-24. 
J 


UNCLASSIFIED 


Approved for public release; distribution is unlimited 


A STUDY OF SECOND AND THIRD ORDER MODEES, FORT tins 
TRACKING SUBSYSTEM OF A RADAR GUIDED MISSILE 


by 


. John Walter Williams 
Captain, СОп:сееаи tartes Пас ше "ОВЕ 
B.S., United States Naval Academy, 1978 


submitted in ракта ртт тшеп ко ЕВС 
requirements for the degree of 


MASTER OF SCIENCE VIN ENSECTRICARhGSEDOCDNIERIDIG 


from the 


NAVAL POSTGRADUATE SCHOOL 
June 1988 


ABSTRACT 


This thesis is a study of missile and target parameters 
used in second and third order modeling of the tracking 
subsystem used in radar guided missiles. Guidance methods 
are analyzed to determine which method is optimum ina 
search for an "ideal" missile. Target parameters which have 
an effect on the missile tracking system are analyzed anda 
target acceleration probability model is discussed. A two 
dimensional third order tracking model is simulated 
utilizing a Kalman filter for target parameter estimation 
and prediction. Linear second and third order tracking 


models are simulated and compared with the third order 


Kalman filter tracker. This thesis concludes that a 
proportional navigation guidance method, with a non linear 
third order tracking Kalman filter, is the better model. 


Benefits of using a non linear third order Kalman filter may 
not overide the cost and complexity of implementation of the 
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qe INTRODUCTION 


Aviation plays an extensive role in current combat 
scenarios. An aircraft, because of its capability to carry 
missiles, is avery formidable weapon platform. Missiles 


provide offensive killing power which change tactics in 


battle scenarios. In order to have the edge in the air to 
air arena, an aircraft must possess the best defensive and 
offensive capabilities. One of the main weapons in the 


aircraft arsenal is the radar guided missile. 

Radar guided missiles are all weather capable, can be 
launched outside of visual range and are less susceptible to 
countermeasures, compared with other missile types. 
Improvement of the missile is a constant necessity to 
Maintain air superiority. 

Improvements in aircraft maneuverability dictate the 
need for missiles to increase performance and capabilities. 
A rule of thumb for design is that the missile must have a 
4:1 acceleration advantage over the target. With modern 
aircraft able to sustain lateral accelerations of ten times 
the force of gravity (G) the missile must be capable of 40G. 

Guidance methods are chosen which optimize the missile 
capabilities to destroy the target. The guidance method 
selected has a large impact on the design of other missile 
subsystems. For any missile to guide to the target, the 
sensor subsystem must track the target. To optimize the 
guidance and tracking of radar guided missiles a predicting 
filter can be used. 

One of the simplest missile guidance techniques 15 а 
second order system which compensates for bearing error and 
bearing rate error. This thesis will look at third order 
models to help optimize the missile sensor subsystem to 


provide better guidance command inputs. A major impetus for 


finding an optimum predicting guidance method is to improve 
missile performance in the final guidance stages. 

There is a region at the end of the missile flight path, 
where the time to intercept is so short that inputs to the 
missile control surfaces will not be effective. If the 
missile has an exact solution of target parameters, it can 
predict the future target position, through the time where 
no inputs will be effective. The missile is flying to the 
projected position of the target, compensating for the time 
delay of control effectiveness. 

Section II will look at some guidance methods and target 
parameters to aid in finding the optimum missile guidance. 


Section III looks at target models and how to implement them 


in the computer simulation. Section IV derives the second 
and third order models. A Kalman Predicting Filter is 
discussed in Section У. Computer simulation and 


implementation of two dimensional third order models with 
the Kalman Filter is given in Section VI. Section VII 
contains conclusions and recommendations. Program listings 


of the computer simulation are in Appendix A and Appendix B. 


five THE IDEAL MISSILE 


Definition of an ideal missile is very difficult. Cost 


or performance functions can be generated to account for 


miss distance, fuel expended, control inputs, flight path 
and numerous other parameters. To the operator the ideal 
missile is one that destroys the target. The designer must 


try to account for a variety of scenarios and targets to 
design the optimum missile. Tradeoffs of performance and 
costs will dictate what the final sub-optimum missile will 
be. 

In an attempt to find an ideal missile insight may be 
gained by looking at the various flight paths and parameters 
for different guidance methods. Three methods to look at 
are pure pursuit, lead pursuit and proportional navigation. 
Pure pursuit would entail a missile always flying directly 
at the target. Lead pursuit would be the case of a missile 
always flying to a point slightly ahead of the target. The 
magnitude of the lead may vary from one time constant, to 
the total missile flight time, ahead of the target. The 
latter would produce a lead collision intercept. 
Proportional navigation utilizes line of sight rate to guide 
the turn of the missile, to zero any further line of sight 
rate. The commanded acceleration to the missile is given by 


the equation: 


АМ = k P Vc (1) 
where AM = missile acceleration 

k = constant of proportionality 

В = antenna angle rate 

Ve = closing velocity 


The constant of proportionality is determined by the 


designer. Observing the effects of the proportionality 


3 


constant on the line of sight rate shown in Figure II-1, any 
constant above 3 is an appropriate value. For k less than 3 
the line of sight rate has its largest slope at the end of 
the intercept. For k greater than 3 the line of sight rate 


will be very small prior to the impact point. 


Line of sight rotational rate 9 
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Figure II-1. Prop Nav Proportionality Constant [Ref. 1] 


Three guidance methods are analyzed to compare flight 
paths, heading changes, line of sight angle and line of 
sight rate. A fourth missile is simulated and plotted, 


identified by "direct". 


The direct path missile is programmed using "a 
posteriori" knowledge of the target flight path. The direct 
path missile goes to the point of impact, in a straight 


line, from the point of launch. 

If the "ideal missile" is defined such that there is no 
missile maneuvering, burns minimal fuel, has the greatest 
launch distance, minimum intercept time and accounts for all 


target maneuver, it will be a direct path missile. 


4 


Two programs were written using Dynamic Simulation 


Language, DSL, to analyze the guidance methods and produce 
plots. The programs are included in Appendix A. 

Three scenarios are used to compare the guidance 
methods: 


- head-on aspect 
- tail aspect 
- beam aspect 
Three target accelerations are used to determine the 
effect of the target maneuver on the parameters. The 
selected target accelerations are: 
ааа с 
WF G'S 
=O. gd 5 
A second order proportional navigation missile with a 
proportional navigation guidance constant of four, where 
the line of sight angle and angle rate is estimated by the 
antenna parameters of angular position and angular rate. 


The governing equations for tracking in bearing are: 


т. 

в | ђак + р, (2) 
© 

ве пра + ñ (3) 
© 

В = -20*B + 100*(LOS-B8) (4) 


antenna position angle 


£ 
ыу 
(D 
H 
(D 
тю 
II 


В - antenna angle rate 
В = antenna angular acceleration 
LOS - actual target bearing 


Pure pursuit and lead pursuit guidance missiles are 
initialized heading directly at the target. Pure pursuit 


guidance maintains heading directly at the target by 


S 


Calculating the heading 


at each step of the discrete 


Simulation. The equation for pure pursuit heading is: 
PPHDG = atan((yt-ym) / (xt-xm) ) (5) 
where yt = current target Y coordinate 
xt = current target X coordinate 
ym = current missile Y coordinate 
xm = current missile X coordinate 


Lead pursuit maintains 


a heading in front of the target, 


with a variable lead, calculated using half the time to go, 


given by: 
LPHDG = atan((yt+tvely*.5ttg)/(xt+tvelx*.5*ttg)) (6) 
where Nt - current target Y coordinate 
XE = current target X coordinate 
tvely = target velocity in Y direction 
tvelx = target velocity in X direction 
tig = time to go in the intercept 


The results of the 


comparisons are given in graph form 


in Figure II-2 through PrYgure s ER. M 


A.  HEAD-ON ASPECT 


Figures II-2 through II-13 are the results of missile 


guidance comparisons for head-on aspect initial condition 


with 0G, 3G and 6G constant target acceleration. The 


missile begins at the 


initial position is at x 


origin of the graph. The target 
= 10000 ft, у = 500 ft. Applied 


lateral target acceleration is away from the missile. 
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Figure II-3 Missile Heading Plot for Head-on 


Aspect No Target Turn 
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Figure II-4 Line of Sight Àngle Plot for Head-on 
Aspect No Target Turn 
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Figure II-5 Line of Sight Rate Plot for Head-on 


Aspect Nc Target Turn 
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Figure II-? Line of Sight Angle Plot for Head-on 
Aspect Constant 3G Target Turn 
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Figure II-10 Position Plot for Head-on 
Aspect Constant 66 Target Turn 
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Figure II-11 Missile Headings Plot for Head-on 
Aspect Constant 6G Target Turn 
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Figure 11-12 Line of Sigh* Angle Plot for Head-on 
Aspect Constant 6G Target Turn 
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Figure II-13 Line of Sight Rate Plot for Head-on 
Aspect Constant 6G Target Turn 
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Figure II-15 Missile Heading Plot Tail 
Aspect Constant No Target Turn 
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Figure II-16 Line of Sight пыте Plot Tail 
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Figure II-17 Line of Sight Rate Plot Tail 
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Figure] 20 Line of Sight Angle Plot Tail 
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Figure I1-21 Line of Sight Rate Plot Tail 
Aspect Constant 3G Target Turn 
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By applying the acceleration away from the missile the 
pilot will lose sight of the missile. This is undesirable, 
but a turn into the missile will help the missile, in the 
early stages, more than a turn away. With a turn into the 
missile, the pilot will also lose sight of the missile, 
during a constant acceleration turn. 

1. OG Target Acceleration 

From the position plot, Figure II-2, it is seen that 
proportional navigation guidance is essentially the same as 
the direct path guidance. Any errors are due to 
initialization of the heading for the proportional 
navigation guidance. The pure pursuit guidance missile and 
the lead pursuit guidance missile fly curvilinear paths to 
targetintercept. The curve for the lead pursuit guidance 
missile is less than the pure pursuit missile due to target 
lead. 

The missile headings graph, Figure II-3, shows the 
relative heading changes involved for each missile guidance. 
After the initialization errors have been corrected, the 
proportional navigation guidance missile parallels the 
direct path missile. The heading changes for the lead 
pursuit are less than for the pure pursuit guidance method. 
Large increases in missile headings at the end of the 
intercept implies large lateral accelerations are required 
for the missile to complete the intercept. 


The line of sight graph, Figure 11-4, shows what 


would be expected for this case. The proportional 
navigation guidance and direct path missile maintain 
constant line of sight, approximately, while the line of 


Sight increases for lead pursuit and pure pursuit guidance 
methods. The large change in line of sight at the end of 
the intercept also correlates to a high lateral acceleration 
required by the missile. 

The line of sight rate graph, Figure II-5, gives 


some insight to the control inputs to the missile guidance 
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subsystem. The values from, Figure II-5, are the slew rated 
for the sensor subsystem. A positive slew rate is seen for 
the prop nav guidance only at the final stages of the 
intercept. Lead pursuit and pure pursuit guidance methods 
have accelerating positive slew rates throughout the 
intercept. The direct path missile has a negative slew 


rate, caused by the missile speed advantage (2500 : 667 
ft/sec). 
2. 3G Target Acceleration 


The position graph, Figure II-6, shows a curvilinear 
path for all three guidance methods. The curvature of the 
target flight path is misleading, because of the axis 
scaling. The target is maintaining a constant acceleration. 
All three guidance methods appear to end up in a tail chase. 
The scaling is misleading again. Target heading change is 
approximately 50 degrees. The proportional navigation 
missile impacts in the beam while lead pursuit will be rear 
quarter and pure pursuit will be a tail chase. 

The missile headings graph, Figure II-7, shows the 
proportional navigation guidance has the lowest heading 
Slope and 1S approximately linear at the end of the 
intercept. Pure pursuit and lead pursuit guidance methods 
have accelerating missile heading slopes requiring higher 
missile acceleration. 

The line of sight graph, Figure II-8, is similar to 
the Figure II-4, proportional navigation guidance method, 
which has low line of sight angles, slightly increasing due 
to target acceleration. Lead pursuit and pure pursuit 
guidance methods have line of sight angles which increase at 
an accelerated rate throughout the intercept. The direct 
missile has decreasing line of sight. The large negative 
LOS for the direct missile at the end of the intercepts is 
caused by the miss distance and heading initialization. 

Line of sight rates, Figure II-9, correlate with the 


line of sight plot, Figure II-8. Line of sight rates are 
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small but show an acceleration at the end of the intercept, 

due to decreasing range. Proportional navigation guidance 

methods are reducing the line of sight angle while lead 

pursuit and pure pursuit increase the line of sight angle. 
3. 6G Target Acceleration 

Comparing the position plot, Figure II-10, with that 
of the 3G case, Figure II-6, similar statements can be made 
about all of the missile paths. Scaling is slightly 
deceiving; the target has made approximately 100° heading 
change. All flight paths are curvilinear with the 
proportional navigation guidance method being the shorter of 
the three methods. 

Figure II-11, shows the heading changes for the 
missiles and the smaller missile maneuvering required for 
the proportional navigation missile. Figure 11-12 апа 
Figure II-13 show larger magnitudes for line of sight angle 
and line of sight rate than the 3G case, but follow the same 
trends. The direct path missile shows a reversal in line of 


sight rate as the target heading change is greater than 909?. 


BEL TAIL ASPECT 


Figure II-14 through Figure II-25 are the results of 
missile guidance comparisons for tail aspect initial 
condition with OG, 3G and 66 constant target acceleration. 
The missile begins at the origin of the graph. The target 
initial position is X=10000 ft. and Y=1000 ft., with an 
initial heading of 090, parallel to the X axis. Applied 
target acceleration is directed into the missile, 
perpendicular to the target heading. 

1. OG Target Acceleration 

The position plot, Figure 11-14, shows the 
proportional navigation guidance missile flies a similar 
path as the direct path missile. The difference in the 


flight paths is due to errors in initialization. The 
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missile heading plot, Figure 11-15, shows that the heading 
for proportional navigation guidance and direct path 
Missiles are parallel after the initialization errors are 
corrected. Pure pursuit and lead pursuit guidance methods 
have continually changing headings with accelerating slopes 
at the final stage of the intercept. 

Line of sight angles for the proportional navigation 
guidance and direct path are approximately constant and 
equal to the initial line of sight angle, giving a constant 
bearing decreasing range trajectory as seen by the target. 
The line of sight angle for pure pursuit and lead pursuit 
guidance decrease, but non linearly, as seen in Figure II-16 
and Figure II-17. 

2. 3G Target Acceleration 

With target acceleration, all three missiles fly a 
curvilinear path. The proportional navigation guidance 
method has the shortest flight path, as seen in Figure II- 
19. Proportional navigation guidance gives a linear heading 
change, as seen in Figure  II-19. Pure pursuit and lead 
pursuit guidance methods give higher heading slopes when the 
target applies lateral acceleration as compared to the OG 
heading plot, Figure II-15. 

Line of sight angle changes are small for propor- 
tional navigation guidance, as shown in Figure II-20. Pure 
pursuit and lead pursuit guidance have decreasing line of 
Sight angles, with corresponding decreasing line of sight 
rates, as seen in Figure II-20 and Figure II-21. Line of 
sight rates increase for proportional navigation guidance, 
as would be expected from the path the missile flies. The 
direct path gives both a nonlinear line of sight angle and 
line of sight rate throughout the intercept. 

3. 6G Target Acceleration 

When the target acceleration is increased, flight 

paths have a larger curvature, as seen in Figure II-22. The 


scaling gives some distortion, the target has gone through 
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approximately 100? of heading change. The missile heading 
changes, as per Figure II-23, are similar to those observed 
for the 6G head-on aspect, Figure II-11. 

The line of sight angle and line of sight rate are 
larger for an increase іп lateral target acceleration, as 
seen by comparing Figure II-24 and Figure II-25 with the 3G 
case, Figure II-20 and Figure II-21. The line of sight and 
line of sight rate for the direct path have a very large 
slope at the final intercept due to effects of decreased 
range. The reversal of line of sight rate for the direct 
path in Figure II-25 is where the target heading change is 
20°. 


C. BEAM ASPECT 


Figures II-26 through II-37 are the result of missile 
guidance comparisons for beam aspect initial conditions with 
OG, 3G and 6G constant target accelerations. The missile 
begins at the origin of the graph. The target initial 
position is X=15000, Y=0. Applied acceleration is directed 
into the missile. 

1. OG Target Acceleration 

As in the two previous cases, with no lateral target 
acceleration, proportional navigation guidance and direct 
Path missiles have similar flight paths, Pure pursuit and 
lead pursuit guidance have curvilinear flight paths, as seen 
in Figure II-26. Heading changes are small for proportional 
Navigation guidance and direct flight path missiles but not 
zero as what might be inferred from Figure II-27, because of 
Scaling. The heading change for pure pursuit and lead 
pursuit guidance is accelerating throughout the flight time 
with the intercept ending in a tail chase. 

Line of sight and line of sight rate, Figure II-28 
and Figure II-29, are similar to the two previous cases, for 


no target acceleration and the analysis is the same. 
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2: 3G Target Acceleration 
For the beam aspect initial condition, when lateral 


acceleration is applied, the difference between flight paths 
for proportional navigation guidance and direct path is 
opposite from the two previous cases for lateral 
acceleration. The proportional navigation guidance missile 
flight path is on the opposite side of the direct path from 
pure pursuit and lead pursuit guidance flight paths, Figure 
II-30, with the opposite curvature. Headings for 
proportional navigation guidance continually decrease while 
pure pursuit and lead pursuit guidance increase. 

Differences between the methods are enhanced by the 
line of sight angle and line of sight rate plots in Figure 
II-32 and Figure 11-33. Proportional navigation guidance 
decreases line of sight while the others have an increasing 
line of sight and appropriate line of sight rate. 

3. 6G Target Acceleration 

Increasing the target lateral acceleration magnifies 
the flight path differences between the guidance methods. 
As has been seen from the previous cases, the larger lateral 
acceleration increases the magnitudes of the values for 
Figure 11-33 through Figure II-37, compared with similar 
graphs from the other cases, but the trends remain the same. 
Proportional navigation guidance parameters have smaller 
changes than pure pursuit and lead pursuit guidance, with 
parameters generally decreasing instead of increasing for 


pure pursuit and lead pursuit guidance. 


D. CONCLUSIONS 


For scenarios with no applied target lateral 


acceleration the proportional navigation missile is the same 


as the direct missile. The pure pursuit and lead pursuit 
missiles finish іп а ‘tail chase' where a missile speed 
advantage is required to complete the intercept. The line 
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EN Sight remains constant for the proportional navigation 
missile but increases with the pursuit missiles. The line 
of sight rate increases with decreased range for the pursuit 
missiles but is zero for the proportional navigation and 
direct missiles. 

When target lateral acceleration is applied, there is a 
deviation between the direct and the proportional navigation 
missiles. Since the target is turning into the missile, the 
line of sight angle decreases at an accelerated rate as 
range decreases. The proportional navigation missile 
accounts for the change of line of sight by turning into the 
target. The pursuit missiles fly a tail chase profile with 
higher line of sight accelerations due to the target turn. 

For the direct missile, when target acceleration is 
applied, the line of sight is not constant and the rate of 
change depends on the applied acceleration. Implementation 
of а direct missile is impossible because the parameters 
used to guide the missile are dependent on the target flight 
path. 

If an optimum missile is to be designed, proportional 
Navigation guidance is the closest to an "ideal" missile. 
The better the proportional navigation missile can 
compensate for the effects of the target acceleration, the 


closer to "ideal" the missile will become. 
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TIT. TARGET IMODEL 


A complete target model for use in computer simulation 
is very involved, time consuming and computer intensive. To 
simplify target simulation the target flight profile is 
based on the fact that the missile sees only the effects of 
the target command inputs and resultant flight path. The 
target model was simplified to include only the flight 
profile desired and not be concerned with the full target 
modeling. A constant speed, constant acceleration target is 
assumed for the simulation. A variable speed, variable 
acceleration target can be added at a later time. The 
missile simulation estimates and predicts target parameters 
of range, range rate, range acceleration, bearing, bearing 
rate, and bearing acceleration. Therefore, for proper 
evaluation of the missile guidance and missile flight 
profiles the target parameters in missile coordinates for 
all these parameters must be computed. 

Complete analysis of target motion is obtained from a 
three dimensional derivation, but insight can be gained from 
two dimensional modeling. Three dimensional flight profiles 
are easily implemented on the computer but graphic display 
of the results are difficult. Two dimensional displays are 
easier to implement and comprehend. A two dimensional 
target model is assumed. 

A target can accelerate at values ranging from negative 
maximum instantaneous acceleration to positive maximum 
instantaneous acceleration, A(max inst). A(max inst) is 
defined as the aerodynamic acceleration given by the maximum 
deflection of control surfaces. A(max inst) is dependent on 
airspeed and air density. High speeds and low altitudes 
produce the highest instantaneous accelerations. Maximum 
sustained acceleration, A(max sust), is defined as the 


aerodynamic acceleration to maintain constant airspeed and 
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Esnstant altitude, at full thrust. A(max sust) can be 
exceeded but must be compensated for by a reduction in 
alrspeed or altitude. 

In three dimensional maneuvering cross coupling exists 
between horizontal and vertical angle and angle rate 
components of target velocity and target acceleration. 
Applied accelerations and velocity changes in one direction 
will affect the parameters, seen Бу а missile, in the two 
other directions. 

Thrust capabilities have a direct correlation to A(max 
sust) and the airspeed of an aircraft. An aircraft with 
higher thrust can maintain a higher speed and compensate for 
drag induced by the applied acceleration. A modern aircraft 
with a relatively high thrust to weight ratio will have a 
very high A(max sust) which is close to A(max inst). 

Airspeed is a key element for maneuverability and 
Survivability. Tactics incorporate optimum techniques for 
Maintaining airspeed or recovering lost airspeed. Pilots 
learn to compensate for limitations of A(max sust) by 
intentionally decreasing altitude and use the effects of 
gravity to maintain airspeed when lateral acceleration is 
applied. Another technique is to apply the lateral 
acceleration required to perform a maneuver then to reduce 
the acceleration, allowing excess thrust to restore the 
airspeed and altitude lost during the maneuver. There is a 
recovery time for the thrust to restore the lost energy, so 


to aid in restoring airspeed, a pilot will normally go to 


zero acceleration, reducing any induced drag, effectively 
increasing the aircraft thrust. This maneuver causes a loss 
in altitude due Lo gravity but improves  airspeed 
restoration. 


The probability of aircraft acceleration is used to 
determine parameters for the target model used in missile 
designs and simulations. Figure ІІІ-1 shows a typical 


acceleration probability graph used in missile design. The 


E 


figure does not account for pilot tendencies nor the 
difference between A(max sust) and A(max inst). The graph 
assumes that а pilot will maneuver primarily at zero 
acceleration, straight and level, or maximum acceleration, 
for a turn, with some probability for any other possible 
acceleration. The design engineer assigns probabilities for 
the impulse functions at zero acceleration and at maximum 
acceleration depending on the type of target aircraft. A 
large bomber may have an A(max sust) half that of a fighter 
aircraft, with less probability of turning than flying 


straight and level. 


PCa? 


-А 
ES acceleration E 


Figure III-1 Probability of Aircraft Acceleration [Ref. 2,3] 


A proper target maneuver model should include some pilot 
tendencies and known tactics. A pilot is not always able to 
move the control surface to a precise location to cause a 
precise acceleration at an optimal time. A pilot will move 
the control surface, judge the acceleration induced then 
move the control surface to achieve a desired acceleration. 
The feel a pilot receives from the ‘'stick' is a prime 


feedback source to allow the pilot to set the desired 


acceleration. The more force the pilot applies, the more 
the control surface moves, and the higher is the 
acceleration. Modern aircraft may use computers to achieve 
the commanded acceleration, reducing any pilot induced 


errors on input. 
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TACTICS FOR MISSILE DEFENSE 


When a missile is fired at the target the battlefield 
scenario changes to an immediate survival situation. If the 
missile is undetected, the acceleration probabilities of 


Figure III-1 may be an adequate target model. If a pilot 


sees the missile, pilot reactions will change the 
probabilities. How the probabilities change may be of 
consequence to the missile guidance. An optimum missile 


design may be able to use pilot tendencies to increase 
missile performance. The overall acceleration probability 
from Figure III-1 is a zero mean function with a variance, 
o2, dependent on the probabilities assigned. Reference 2 
discusses obtaining the parameters for the target accel- 
eration probability model of Figure III-1. 

When the pilot imposes amissile defense, the overall 
acceleration may not be zero mean, nor maintain the same 
variance. To account for the changes in acceleration 
probability a function similar to Figure III-2 might be 
used. This model accounts for some variance to the 
acceleration which the pilot is trying to achieve centered 
around A(max sust). If the pilot is trying to achieve A(max 
inst), it 15 assumed he will be decreasing airspeed and 
reducing actual acceleration until the applied acceleration 
1s decreased to A(max sust) or below. A smart pilot will 
either fly at a maximum acceleration or at zero 
acceleration, increasing the aircraft maneuverability. 


Last ditch maneuvers are performed at A(max inst) to 


avoid the missile, neglecting any adverse effects of 
applying acceleration, in order to increase survivability. 
If the last ditch maneuver is performed too soon, 


acceleration is decreased, due to loss of airspeed, negating 
the effectiveness of the maneuver. Further studies can be 
made correlating the use of A(max inst) versus A(max sust) 


for missile defense tactics. 


Sd 


P(a) 


0 max max 
sust inst 
ACCELERATION 
Figure III-2 Probability of Acceleration 

If pilot reaction is taken as Gaussian when trying to 
achieve a desired acceleration, the overall acceleration 
probability will be Gaussian with а non-zero mean anda 
variance dependent on the combination of the two Gaussian 
terms. The probability assigned to each term determines the 
mean and variance. 

Missile defense includes placing the missile on the 
beam, to utilize the largest acceleration vector, with 
maneuvers made out of phase, out of plane with the missile. 
The largest acceleration component comes from the elevator, 
perpendicular to the wings. By placing the wings parallel 
with the plane of the missile, the largest acceleration 
component is used to create the largest missile corrections, 
perpendicular to the plane of the missile. The plane of the 
missile is defined by three points: target position, 


missile position and the projected impact point. 


A graph of a target acceleration, while performing 
missile evasion, might look like Figure III-3. The pilot 
commands A(max sust) or A(max inst) for a short time, then 
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reduces the acceleration to zero to regain lost energy, 
before applying the acceleration again. This process may be 
repeated 2,3 or more times during the missile flight time. 
The graph attempts to incorporate transient response induced 
by system time delays, transient responses of the control 


surfaces and pilot tendencies for maneuvering and control of 
the target. 


ACCELERATION (G) 





TIME (sec) 


Figure III-3 Target Acceleration 
The resultant average acceleration is non-zero, with a 
non-zero variance. The mean and variance of Figure III-3 
can be estimated by the parameters assumed in Figure III-2. 


Figure ІІІ-1 and Figure III-2 show total aircraft 


acceleration. The parameters as seen by the missile, in 

antenna coordinates, will vary depending on the three 
dimensional maneuver employed. The missile tracking 
subsystem must be designed to handle the maximum 


acceleration possible for each of the orthogonal components 
of the reference frame. 

Proportional navigation missiles compensate for the 
applied target acceleration by decreasing the line of sight 
rate induced by the change of target velocities. With a 
good target model and detection of maneuvering effects, the 


missile guidance can predict target motion and position. 
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TV. SYSTEM MODEL 


From the section on the ideal missile, it was 
ascertained that target parameters are not constant if 
lateral acceleration is applied. This thesis will attempt 


to incorporate as much sensor information as is available in 


defining the system models for an optimum missile. Current 
radars allow the meaSurement of range, range rate and off 
boresight bearing error. Measurements taken by the radar 

are referenced to a radar axis system. In order for the 


missile to use the information supplied by the radar, a 


common reference frame must be established. 


A. COORDINATE SYSTEM 


Each entity in the missile-target intercept problem has 
its own coordinate system. The overall geometry as seen 
from an "eye in the sky" would view it in space coordinates. 
An observer on the ground would view it in earth 
coordinates. The launching aircraft, missile and target 
aircraft will view it іп ап individual coordinate system, 
referenced to that specific platform. Trying to equate each 
coordinate system is not an easy task but one which is done. 

By use of Euler angles any reference coordinate system can 
be related with Earth coordinates. By use of a directional 
cosine matrix transformation any reference coordinate system 
can be transformed to another reference coordinate system 
[Ref. 3]. 

The missile is concerned with flying to a point in space 
that will hopefully be occupied by the target at the 
completion of the intercept. The object is to guide the 
missile to the proper point where the target will be. The 
missile is concerned with its coordinate system and not 


that of the target. But on the missile itself there are 
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various coordinate system reference points. Each sensor has 
its own location on the missile and where it is mounted is 
its reference point. Any moving sensor, like the antenna, 
will have its special reference coordinate system. Missile 
parameters are normally referenced to the missile body frame 
of reference while target parameters are referenced to the 
antenna frame of reference. While very complicated, the 
frames of references can be transformed and equated. [Ref. 3] 

To simplify Simulation and evaluation of desired 
parameters, an inertial frame of reference will be used 
which is centered at the radar antenna location. This 
simplification will aid in better evaluation of the effects 
of the target parameters andthe missile guidance without 
encumbrance of transformation errors. Although the 
Simplification assumes ideal missile parameters, time delays 
can be incorporated later to account for first order 


modeling of the missile. 


B. EQUATIONS OF MOTION 


In cartesian coordinates missile and target motion is 


described by the standard motion equation: 


T T 
X(T) = Xo + | x(t) at + || iit) dt at (7) 
O O 
The equation is based on a fixed reference point. The 


orthogonal directions (Y and Z) will have the same equation. 
The antenna frame of reference uses polar coordinates 


which have the equations: 


T. TS ~. 

r(t) = ro + | r(t) dt + Ц r(t) dt dt (8) 
о о 
те тт. 

e(t) = e. + | e(t) dt + || e(t) dt dt (9) 
о о 
"TM E 8 

$(t) = Ф. + | @(t) dt + || $t) dat dt (10) 
о о 
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where r(t) » the radial component of motion 
@(t) = horizontal component of motion 
@(t) = vertical component of motion 


When the reference point is not fixed, extra terms and 
cross coupling are introduced into equation dynamics. The 
coriolis equation accounts for the moving reference point. 
Depicted in Figure IV-1 a change in the vector К 15 
accounted for by both changes in the magnitude and the 
rotation effects by the moving reference point. Using the 
terms as defined in Figure IV-1 we can obtain the necessary 
equations to find r(t), @(t) and @(t). For simplicity, only 
the derivation for r and © are shown with r and Ф 
relationships being a duality of derivation of r and 8. The 


Simulation of Section VI will be two dimensional. 





Figure IV-1 Rotating Vector Diagram 
The coriolis equation to relate the time rate of change 
of the R vector to the rate of change in the r direction and 


the angular rate of motion is: 


К = r + wxR (11) 
where R = the directional vector 

r = the magnitude of the directional vector 

w = the angular rate of motion 
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The total change of the vector R is the sum of the change in 
the magnitude of R due to changes along the original vector 
R given by r and the angular rotation due to the moving 
coordinate frame, given by wxR. 

Utilizing the general rule for differentiating a vector! 


an expression is obtained for the acceleration of the R 


vector. 
R = r + wxr + wxr + wxR (12) 
R = г + wer + ихг + ихг + ихихг (13) 
R = r + wxr + 2(wxr) + ихихг (14) 
where R = directional vector 


r the acceleration of the magnitude of R 
wxr and wxr are perpendicular to the R vector 


WXWXr is centrifugal acceleration 


This equation gives the cross correlation of range and 
angle to implement in a second order model. Applying the 
rule of differentiating a vector again will yield the 


equations for a third order model. 


d(R ) = r + w xr + Wxr + 2 (wxr) + 2(wxr ) * WXWXY 
at : 5 ҚҰ! 
+ ихихг + wxwxr + wxR (15) 
A(R) = r +w xr + wxr + 2 (wxr) ~ 2(wxr ) + ихихг 
dt * е v е е е 
+ WXWXI + WXWXYI + WXI + ихихг + мх2 (мхг) + MEI 
R = r + 3 (wxr) + 3 (ихг ) + W Xr + мхмхг 
+ 2 (wxwxr) + 3 (wxwxr) + ихихихг (17) 


where R is the acceleration jerk of vector R 


1 rule for differentiating a vector 
= = 2A + WXA 
9 
the total derivative is the sum of the time rate of change 
of the vector and the rotation of the vector. 
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r is the magnitude of R 
Wxr is perpendicular to R 
wxwxr is in the negative direction of R 


WXWXWXr is perpendicular to R 
Es SECOND ORDER MODEL 


Beginning with equations 8-10, а second order state 


Space system can be set up which would have the form: 


: T 
r= Yo + E r dt + || r dt (18) 
O O 


т. 607 
Ө = Ө, + | e dt + || e at (19) 

O O 

r 0 1 0 0 Y 

r о о о 0 r 

MP (20) 

e ОАО О Ө 

е оо оро е 


The range portion of the second order system may be 
accomplished totally by the radar, Since no other subsystems 
require the information. The radar receiver is designed to 
track the target in the radial direction without the need 
for an additional filter. 


The more difficult state equations to implement in the 


missile are the angular directions. The second order, time 
invariant, constant velocity, zero acceleration, state 
feedback model makes the tracking much easier. The 


continuous system model can be given ina time derivative 


form as: 
X + ki-X + kz -X = 0 (21) 


If the term given in the equation as X is actually the error 
of the angular position then 


X - (line of sight angle - antenna position angle) 
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where X can be directly measured by the antenna. The values 
of the k's in the time derivative equation depend on the 
designer and the response desired. For the simulations of 
the proportional navigation missiles of Section II and 
Section VI, ki = 20 and kz = 100 were used. These constants 
give a response time constant of 0.1 sec. 

The use of the coriolis equation to derive the second 
order model gives a time varying solution. Implementing 
time varying equations are difficult and often avoided by 
using a time invariant system and state feedback to cancel 
errors. The time variant space state model derived from 


equation 14 is shown below: 
В = г + wxr + 2(wxr) + wxwxr 122) 


Separating into orthogonal components of radial and 


transverse with scalar multiplication: 


II 


ав r + миг (23) 


ar = wr + 2(wr) (24) 


Rearranging into equations to implement into a system: 


r= -@2r + ar (25) 
6 - -20r + ar (26) 
IÇ r 


The state space model is difficult to represent unless 


divided into two channels with cross coupling given in the A 


matrix. 

т 0 1 r 0 ак 
к = : . + (27) 
2 -82 0 I 

9111 e aT 

= Е , + (28) 
0 -2r e 
r 
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D. THIRD ORDER MODEL 
The time invariant model, derived from equations similar 


to those deriving the time invariant second order model 


(equations 18 and 19), in state space form is given by: 
r 0 1 о о 0 O r 
r` 0 0 1 0 0 O r 
г [= [о о o о о о r (29) 
8 0 0 0 0 1 O e 
е 0970 0 00 4 е 
е ооо ООО е 


Ranging may be accomplished by the radar receiver, as in 
the second order model, for similar reasons. Tracking 
angular positioning requires knowledge of the target angular 
acceleration аз well as angular velocity. À filter is 
normally used to maintain a track of the target angular 
parameters. Common filters are a-pB, Weiner, and Kalman. As 
compared in Reference 4, the Kalman is the best filter 
suited for air to air missiles, but also the most costly to 
implement. For the time invariant third order model a 
Simple constant gain Kalman Filter can be used. The Kalman 
Filter will be discussed in the next section. 

The tine variant third order model is obtained from the 
second derivative of the coriolis equation derived in the 


previous section. 


`R Y + 3(wxr) + 3(wxr ) + w xr + WXWXY 


+ 2(wxwxr) + 3(wxwxr) + WXWXWXI (30) 

The R term is the change in the acceleration of the 
vector or a "jerk" term, a simple comprehension is the rate 
at which the pilot applies the commanded acceleration. 
Separating the equation into radial and tangential terms, 


the two orthogonal scalar equations are: 
ae = г - 3w2r - 3wwr (31) 
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ar = 3Wr + 3wr + w r + w?r (32) 
Converting the equatlon into primary coordinate axis, 


the equations obtained are: 


г = ар + 3з@2г + 39(0 )г (33) 


6 = ar - 3(@ )r - зе(г) - ез (34) 
Y r 


The disassociated space state model looks like: 


0 1 г ак 
- 0 0 1 r + 0 (35) 
зе(е ) 382 r 
е 0 1 0 е ат 
ө |= 0 0 1 е + (36) 
6 0 -3r  -3r e 
r Е 


It is readily observed that a very high cross coupling 
of the radial and transverse components exists. The range, 


range rate and range acceleration are required to adequately 


compute the angular acceleration. The angular velocity and 
acceleration is required in computing the range 
acceleration. All of these quantities are time varying 


requiring a time varying filter to implement this model. 

The cubic term of angle rate inthe angle channel is 
insignificant compared to the other terms and is neglected. 
The second order model uses the simplifying assumption of 
constant velocity and constant acceleration. For the full 
third order model, no simplifying assumptions will be made. 
This third order model should account for all of the cross 
coupling between the bearing and angle channels. 

A Kalman Filter can be employed to track the target in 
both range and bearing to implement the time variant third 
order model. The Kalman Filter will be discussed in the 


next section. 
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V. KALMAN FILTER 


Given a system model, where the plant can be modeled by 
a set of first order differential equations and the output 
can be measured, a set of state equations can be defined 


similar to: 


X = АХ + ВИ +Н (37) 
Y = C X + V (38) 
where X is the state vector 


Y is the system output vector 

A,B,C and D are matrices 

U is the system input 

W is plant disturbances 

V is measurement noise 
The system can be modeled in discrete time as: 
X(k+1) = @ X(k) * T U(k) + W(k) ( 
Ки HOOK) SENE (40) 
À Kalman Filter is the best filter to track the output 

of а discrete system [Ref. 3]. The Kalman Filter equations 


are given as: 


X (w |x) = отказ) + G(x). | Yo) = "TE || (41) 
x p) - $-X( p) + rU (X) (42) 
YO: p) - нохо р) (43) 
where KO [i = the state estimate at time k given 


К information through time k. 
KA(**? |e) = the state estimate ас time ki an 
Е information through time k. 

УК еј = реноме Ес станса ел 2 
given information through time k. 


G(k) = the filter gain at time k. 
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For linear, time invariant systems, the d$ and r matrices 
are easy to calculate and follow dar oct ly from the state 
Space model, where Ф = е апа Г = Е dt. For non linear 
Systems, an extended Kalman filter can be used. For the 
extended Kalman filter, the Ф and T matrices are linearized 
about the projected operating point. One method of 
estimating the linearization is to take the partial 


derivatives of the non linear state space matrices: 








Ф = А (44) 
S 
X = Xo 
U = Uo 
и = 
Г = В (45) 
U 
X = Ko 
U = Uo 
W = 0 
The gain matrix G(k) will vary with the parameters of 
the filter. G(k) is the weighting factor for the system 
error. The solution to the filter gain G(k) requires the 
solution of Riccati equations: 
EY 
С(К) < РК i) Hr [8 Bei) H * ROO] (46) 
Р(К |к-1) QPPL фт + ó Q от (47) 
Ek) - Р(Е |к-1) — С{(К) в ЕШ БЕТ!) (48) 
where G(k) + Kalman Filter gain at time k 
Pet eee OvVardance Of predicted estimate 
R(k) » measurement covariance matrix, EÍVV'!| 
Q(k) = maneuver covariance matrix, E[{UU™ } 
НОЕ Co variance of faltered estimate 
6 = maneuvering weighting matrix 


A constant gain matrix can also be used in the Kalman 


filter. Instead solving the full Riccati equations for each 
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change of variables, a constant value is used throughout the 
problem. À constant gain matrix will simplify 
implementation of the Kalman filter. 

One implementation of the third order model, as 
discussed in the previous section, is to model the system as 


linear, time invariant, given by the space state model: 


(49) 


QD. QD. 0) • қ.» қ.» қ. 
| 
о о о о o o 
OPO (ж) оо 
О О О СО ко 
о о о о о o 
ое ео 
о е o o o o 
Ф. QD: D Me He Fs 


The Kalman Filter equations for this third order model 


аге. 
X( j) Gs ü j о ре м (50) 
X(k*i |k) = &:X(C p) (51) 
Уат |к) = H-X(%*1 у) (52) 
where X = E and H = 1 0 070 00 
| 010000 
r 001000 
а 000100 
r 000010 
000001 
е 
ë 
5 


Using a Kalman Filter on this third order model is very 
simple and requires few on line calculations. The gain 
matrix can be considered either constant or time varying. 
If time varying gains are used, they can be computed off- 
line and stored in memory. The Filter then utilizes the 


precomputed gain schedule and can select a gain depending on 
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the accuracy of the filter at that time. If a maneuver 
causes the filter to loose accuracy then a higher gain term 
can be utilized. If constant gains are used then they must 
be high enough to compensate for any maneuver the target 
might make. A high gain matrix will make the missile more 
responsive to any unwanted noise terms in the system since 
the missile cannot distinquish a noise input from a target 
maneuver. 

As discussed previously, a Kalman Filter is not required 
for the range channel. The radar can meaSure range and 
range rate directly. Since the actual values of the range 
channel are not used by any other elements of the guidance 
Subsystem, the radar is able to maintain its own tracking of 
the target in the center of the range gate, which has no 
consequences on the rest of the missile guidance. Some 
noise information can be gained when estimating the range 
channel with a Kalman Filter. A Kalman Filter is used for 
the range channel in the simulation of Section VI for 
completeness of simulation and practical experience. 

A second implementation of the third order model is by 
uSing the equations obtained through the coriolis equations. 


The disassociated state space model is given as: 


0 о г o | ar 
r = 0 ош 1 4 0 (53) 
L 38(8 ) 382 0 r 
e 1 0 е 0 at 
e | = 0 1 8 + 0 (54) 
e 0 ES E в 1 


Two possible ways to implement the Kalman filter are to 
create an extended Kalman filter by linearizing the Ф 
matrix, reducing the time dependence and cross coupling of 


the range and bearing channel, or keep the cross coupling 
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components and have the Kalman Filter maintain the values of 
the time varying Ф. If @ is linearized, some of the cross 
coupling and time dependence lost by linearization can be 
compensated for by the maneuver covariance matrix Q. Given 
the target acceleration probability model, Figure III-1, the 
Q matrix can be calculated, as derived in Reference 3, as 
time varying and relates the cross coupling of the bearing 


and range channel as: 


ОК = O O 0 апа QS = 0 0 (55) 
оо о 0 0 0 
O 0% О O0 0% 
r 
where O2m = acceleration variance 


As discussed in the reference the Q(3,3) element can be 
increased to make the missile gain matrix put more weight on 
any target acceleration elements. 

ТЕ the Ф matrix is maintained as time varying and 
nonlinear then the Q matrix can be constant. The constant Q 


Matrix can be calculated as: 


Qu cc dep T (56) 
where К = matrix dain 
I = identity matrix 


Since the reference deals extensively with the time 
invariant Ф апа the time varying Q matrix, this thesis will 
deal with the time varying Ф and constant Q. 

If Figure III-4 is used to define the maneuver 
probability then a non-zero mean is established. The Q 
matrix maintains the same properties just discussed with a 
different calculation for of The non-zero mean can be 
implemented by increasing Q(3,3), increasing the weighting 


matrix 5, or by not assuming the U(k) term is zero. 
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falman filter equations for this third order model are: 


~ 


E jk) 


^ 


EXE! [y ) 
Yoigo 
Р(К |к-1) 
С (к) 


В к) 


ж 


XC i-i) € GG0-[ YOO ура-а) | 


Ф-Х(К |k) + T-U(k) 
H ENS 
оо о. 0.0: 


Р (к к-а) ент | НАР (к а-а) ент + КОК) | 


Р (к |к- 1) = КЕ) 
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(57) 


(58) 


(59) 


(60) 


(61) 


(62) 


VI. SIMULATION 


To aid in the efforts of simulation, the Dynamis 
Simulation Language (DSL) was used to integrate the 
equations of motion for the target and missiles, as well as 
the antenna angle channel. Two programs were written and 
are listed in Appendix В. The first program is the time 
varying third order model. The second program is the time 


invariant third order model and the second order model. 

The output of the simulation is a set of graphs to 
compare the three missile models and their effectiveness in 
tracking the target. 

The Kalman Filter is implemented in a Fortran subroutine 
at the end of the DSL main program. The basic filter 


equations used were described in the previous section. 


A. ASSUMPTIONS 


The following assumptions are made to simplify the 
implementation of the Kalman Filter and determine the 
effects of the time varying third order model. 

- "Ideal" missile autopilot. 


- Inertial cartesian reference frame for angle 
measurements. 


- Final portion of intercept only. 


- Cross coupled effects of missile motion on antenna 
stabilization system disregarded. 


- Missile initialized to collision course. 


- Missile constant speed of 1500 kts (Mach 1.5) or 
2500 ft/sec. 


- Target constant speed of 400 kts (Mach .75) or 667 
ft/sec. 


- Target lateral acceleration applied perpendicular 
to target velocity vector. 


- Angle of Attack not accounted for in velocity 
vector calculations. 


- Missile located at the center of the cartesian 
reference frame. 
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тате оса с O the right of the origin of the 
reference frame. 


- No noise. 


B. INITIALIZATION 

The user is asked at the beginning of the program to 
establish the geometry by giving the target initial 
position, heading, speed and acceleration. The target 
heading is oriented relative to a vertical line, parallel to 
the Y axis, defining the North or 000 heading. The initial 
missile headings are calculated for constant velocity, zero 
acceleration collision course with a time to go estimate of 
range/missile velocity. Antenna parameters are initialized 


to initial line of sight and zero angular rate. 


C. SECOND ORDER MODEL 

The second order model 15 implemented using a 
proportional navigation constant of four and two s-plane 
poles at $=-10. This gives the equation for angular 


acceleration as: 


Ë = -20 KB + 100 (LOS - В) (ean 
where B = the antenna angular acceleration 

B - the antenna angle rate 

В = the antenna angle position 


LOS = the actual angle to the target 


D. THIRD ORDER MODEL 

The Kalman filter is used to implement the third order 
model. The DSL main program calls the Kalman Filter 
subroutine at a sampling time of h=0.01 seconds. The Kalman 
Filter is executed, then control is passed back to the DSL 
program. A proportional navigation constant of four is used 


as discussed in Section II. 


(1 


1. Time Invariant, Constant Phi Model 
The discrete time invariant third order model 


divided into two Kalman filters for range and bearing is: 


RNG(**! |u) = RPHI-RNG(X |x) (64) 
RNG (X |k) = RNG(*|k-1) + GR(x)+[ DELR | (65) 
S(k*1 |k) = SPHI-S(* jk) (66) 
LO # 509 |к-а) + 65(%). | SDEL | (67) 
where 
RPHI - SPHI = [| 1 .01 .0005 

0 1 

DELR - 


ВВ. РЕ а 
ROM SARO EI 


SDEL 


LOS MSIE NL 
LOSD-=FS ORPI 


Using Matlab functions of Aker and Place, with 


eigenvalues of 0.5,0.5 and 0.5 constant gain matrices were 


obtained for range and bearing, given by GR and сб, 
respectively. 
GR z 0.5 Ооо ов 
0.0025 1.0 
ДОБ 24.9 
GS z QE 
12 5 
1250. 
2. Time Variant, Variable Phi Model 


The discrete, time variant third order model divided 


into two Kalman filters for range and bearing is: 
RNG (k |k) = RNG(* |k-1) + GR(k)-[ DELR | (68) 


RNG (K*! |x) = ВРНТ-ВМС (к [к) (69) 
72 


РК (К |к-1) = ВРНТ-РК (К |к-:)-ВРН1Т + ОК (70) 


-1 
GR(k) - PR(* [x-1) -HRT 7 |BR- PROF [x - 1) «HRT +RMCOV (k] | 
(7%) 
ERU Ik) ее ВНК РЕ: (72) 
S(k |k) = s(k|=-:) + GSQ)-[ SDEL | (73) 
(кат |) = SPHI-S(* |x) (74) 
PS(*|x-1) = SPHI-PS(k |x-1)*SPHIT + QS (75) 
GS(k) = PS [x-1) “HST + [HS+PS(* [x-1) +HST+RSCOV(K) | 
(76) 
PS(*|x) » PS(*|x-1) - GS(X)-HS-PS(* [x-1) (77) 


The time variant model uses two different Ф 
matrices, one for range and the other for bearing. Two 
other matrices must be specified for each filter, the 


Initial error covariance matrix andthe target maneuvering 


covariance matrix. The two € matrices are: 
RPHI = 1 о . 00005 (78) 
.00005%А 1+.00005*B о 
01*A .01*B 1+.00005*B 
SPHI = 16 eu .00005 
0 1*.00005*C 201+. 00C005*p 
0 .01*C+.00005*D 1+.01*D+.00005* (C+D) 
(79) 
where ди- -за (ө) 
В = -зе: 
С = -3r /r 
D = -3r/ r 


The initial error covariance matrices are given as: 

PR (° |o) = 500 0 0 
0 500 0 

0 0 


500 


а 


PS(? |o) = le4 0 0 
0 1е4 0 
0 0 1е4 


The maneuvering covariance matrix can either be time 
varying ог Constant as discussed previously. The 
maneuvering covariance matrix accounts for the capabilities 
of the target as discussed in Section IV. Reference 3 gives 
the derivation for the time varying matrix for a constant € 
matrix. The @ matrix also accounts for any time variance of 
the target model so the maneuvering covariance matrix can be 
constant. A constant matrix is assumed since d is time 
varying. The maneuvering covariance matrix (QR AND QS) are 


given as: 


ОК + 500 0 0 
0 500 0 
0 0 500 
QS = 202 0 0 
0 OT 0 
0 0 01 
The simulation calculates the range model then the 
bearing model. The results of each filter are used in the 


other filter to calculate the ọ values. 
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2. RESULTS 
Similar simulations of the ideal missile cases, Section 
II, were used to evaluate the missiles. The simulation 
consists of head-on, tail and beam aspects with 06, 3G and 
6G target acceleration. The results of the computer simu- 
lation is shown in graph form in Figure  VI-1 through VI-31. 
The computer program listings are contained in Appendix B. 
1. Gains 

Gain comparison plots are given in Figure VI-i 
through Figure VI-4. The resultant gains from the time 
variant, varying Phi third order missile are the same for 
each scenario. 

ВЕСЕ), the varying Phi model 
weights the range error by .5 and the range rate error by 0. 
The constant Phi model uses weights of wand 20125; 
approximately the same, Figure VI-1. 

ОЕ сс КВ}, Lor varying Phi, are 
.001 and .6 while those of the constant Phi are .0025 and 
o Figure  VI-2. Little emphasis is placed on the range 
error, because the radar is measuring range rate, witha 
higher weighting factor. 

In predicting RDD(**! |x), small gains are calculated 
by the varying Phi while the constant Phi model places a 
high emphasis on range rate error. The noise of the system 
will be noticed more in the prediction of КОВАНИ и) than 
the other parameters, because of the higher weighting 
factor. 

Bearing channel gains give unusual curves for the 
varying Phi model, Figure  VI-4. There is little weight 
placed on non-observed parameters as in SG2 and SG3, during 


most of the intercept, except in the initial stage and final 


stage. The gains are highest during the critical stages of 
the intercept.  SG1 is a constant 1.0 giving equal weighting 
to the current estimate and the error. The constant Phi 


model has higher gains giving more weight to any errors. 
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Figure VI-3 Range Covariance Gains RG31 and RG32 
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Figure VI-4 Bearing Covariance Gains SG1, SG2 and SG3 
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Figure VI-25 Missile Commanded Acceleration Plot for beam 
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2.  Head-on Aspect 

Figures VI-5 through VI-16 are the results of 
Missile guidance comparisons for head- on aspect initial 
condition with 0G, 3G and 6G constant target acceleration. 
The missile begins at the origin of the graph. The target 
initial position is at x = 10000 ft, y = 1000 ft. Applied 
lateral target acceleration 1S away from the missile. 

a. OG Target Acceleration 

Flight paths for the three models are shown in 
Figure VI-5. All three paths appear to be the same, within 
the accuracy of the plotter. The line of sight angle, 
Figure VI-6 remains relatively constant for all three 
models, at the initial value, giving а constant bearing 
decreasing range, no manuever intercept. 

Commanded Missile Acceleration, Figure  VI-7, 
shows that the second order model pulls more lateral 
acceleration than third order models. Higher gain terms 
make the constant Phi third order model erratic when 
compensating for initialization error. 

b. 3G Target Acceleration 

In the position plot, Figure VI-8, the second 
order model begins to lag the third order model. The 
lagging means the missile is slower to compensate for line 
of sight rates. This is further illustrated by Figure VI-9, 
the line of sight angle plot. The change in line of sight 
1s greater for the second order model. 

Commanded acceleration, Figure VI-10 shows a 
larger increase for the second order model. The second 
order model requires approximately 7.5G to intercept a 3G 


target while the third order model require 4.66G. 
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c. 6G Target Acceleration 

Increased target acceleration increases the lag 
of the second order model, Figure VI-11. Magnitude of line 
of sight is larger for the second order model, Figure VI-12, 
but the missile is compensating for the errors. The two 
third order models are fairly close. The difference is only 
fractions of radians. 

The constant Phi third order model has less 
initial oscillations in commanded acceleration when the 
target acceleration increases. The second order model 
requires approximately 11G for a 6G target, 7.5G is required 
for the third order models. 

Э Tail ASpect 

Figures VI-13 through Figure VI-21 are the results 
of missile guidance comparisons for tail aspect initial 
condition with 0G, 3G and 6G constant target acceleratior. 
The missile begins at the origin of the graph. The target 
initial position is at x = 15000 ft^ и Е. Applied 
lateral target acceleration is into the missile. 

a. OG Target Acceleration 

As shown in Figure "VI-1353, there is Jitt@e 
difference in flight paths for the three missiles in the no 
target acceleration case. Line of sight angles remain 
constant, Figure VI-14, throughout the intercept 

There is some slight difference in commanded 
accelerations, Figure VI-15, with the second order model 
being the smaller of the three models. 

b. 3G Target Acceleration 

As the target applies acceleration, the lag of 
the second order missile gives a shorter flight path, than 
the third order missiles, Figure VI-16. Line of sight 
angles are smaller for the second order model, Figure VI-17. 

Acceleration required for the second order model 
is 3.4G and 5.3G for the third order models, Figure VI-18. 


By lagging the other missiles, the second order model allows 
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the target to complete part of the intercept, which allows 
the missile to pull less acceleration. 
C. 6G Target Acceleration 
With higher target acceleration, the flight 
paths, Figure VI-19 have approximately the same differences 
as the 3G case. The second order model maintains a better 
flight path throughout the intercept. The slope of the line 
of sight curves, Figure VI-20, are higher with the second 
order model maintaining a smaller angle difference. 
Commanded acceleration is much smaller for the 
second order model than the third order models, Figure VI- 
21. To intercept the 6G target, the second order model 
requires 6.2G and the third order models require 8.1G. 
4. Beam Aspect 
Figures VI-22 through VI-30 are the results of 
missile guidance comparisons for head- on aspect initial 
condition with OG, 3G and 6G constant target acceleration. 
The missile begins at the origin of the graph. The target 
initial position is at x = 15000 ft, y = O ft. Applied 
lateral target acceleration is away from the missile. 
a. OG Target Acceleration 
Figure VI-22 through VI-24, show the three 
missiles are practically the same for a no target 
acceleration, beam ,aspect intercept. А11 three missiles 
maintain a constant bearing decreasing range, small 
acceleration intercept. 
b. 3G Target Acceleration 
Only slight differences are noticed when target 
acceleration is applied. Flight paths, Figure VI-25, shows 
very little deviation. The line of sight angle difference 
of .01 rad between the models is approximately .575 . 
Acceleration is higher for the second order missile to allow 
it to fly the same path as the third order models, Figure 
MI- 27. 
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с. 6G Target Acceleration 


Flight paths have a pronounced difference with a 


higher target acceleration. Pigure VI 22: Line of sight 
angles, Figure VI-29, have larger magnitudes and higher 
slopes. Commanded acceleration increases dramatically, now 


25G is required of the second order model and 9.3G for the 


third order models. 
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Vil. CONCEUS TONS 


In analyzing second and third order missile tracking and 


guidance subsystems, the following conclusions are made: 


dare. 


Proportional navigation guidance is the optimum 
method for missiles, given current design 
tradeoffs. 


Target modelin is very difficult and requires 
the analysis or many factors. Acceleration 
probabilities make modeling easier, but the 
proper acceleration model must be chosen. 


Cross coupling between coordinate reference axis 
components does exist and gives errors if not 
accounted for in the system model. 


Kalman filters are the best predictors for 
alrborne missiles, if one is required. 


Complete time varying third order models give 
better results than approximated linear, time 
зимагдапе third order models. 


Only small differences are noticed in parameter 
values between second and third order models. 
Higher accelerations are required for the second 
order model. 


Second order missiles are better than third order 
missiles in tail aspect, constant acceleration 
intercepts. 


Implementation of a Kalman filter requires | 
considerable amounts of computer resources, with 
limited time to complete the calculations. 


Some parameter terms are of the approximate order 
as system noise or non significant calculations. 


Some recommendations for future study and consideration 


A study of miss distance analysis for second and 
thard order models. 


Analysis of the effects of the Q and P matrix 
Snrcralization. 


Ànalyze the target acceleration probability model 
to find optimum values to assign to the 
probability model. 


Determine missile, cross couple effect of heading 
changes and autopilot torques on the sensor 
subsystem. 


Add noise to the system to determine the effects 
of the high gain terms on a noisy system. 
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BEPBENDIX Агире SS MBE PROGRAM LISTING 


С MISSILE PROGRAM FOR FLIGHT PATH COMPARISON IN THE 
E ео, THIS PROGRAM IS ВЕЗТОМЕП ТО FOLLOW THREE IDEAL 
E ВЕ oe oe OS UNG DUPE ERE NeGUAOANCE TECHNIQUES TO TARGET 
e INTERCEPT. THE GUIDANCE TECHNIQUES ARE: PROPORTIONAL 
e IO U ICGATION PURE PURSUIT AND LEAD PURSUIT. PARAMETERS 
С ARE ALSO OBTAINED FOR A FOURTH MISSILE USING A DIRECT 
e FLIGHT PATH. 
INITIAL 
CONST 
Ы СНЕ ee DE ADEITCHR 4, E17 9.13192, RM1- 20000 
NN=0 
MISSXO=0.0 
MISSYO=0.0 
DN 2500.0 
БОБ 12.5) VT,AT,THDC,TGTXO,TGTYO 
5 HORUM ТЕРІ 25, Е5,.1/2Х,Е6:1,2(2Х7Ғ10.2)) 
TGTV=VT*K2F 
TGTA=-AT*G 
THDG-THDG*DZR 
С 
Е TGTVYO=TGTV*COS (THDG) 


e ТЛ ОТ АТТОО  ОКАОТКЕСТ INTERCEPT MISSILE 
DXT=MISSXO 
DY =o 10 


C 
E ПЛОВЕ TTONOB PROP NAV MISSILE CONSTANT VELOCITY, 
meno ACCEL 
LO O TpOT ом се ХЕ + °(TCGCTYO — MISSYO)**2 )**,5 
TTGO= RO/VM 
ESSE ATAN TCOeOTYO hISSYO,TGTX0=MISSXO) 
BHBG-ATANZ2UCIGTYOTTGTVYO*TTGO-MISSYO,- 
EI XOTIGIVXO*TTGO-MISSXO) 
КОО ТОТУ*^СОС(ТНОС)›/ (СО5${105) % КО) 


E Е LOS 

C а ФАТТОН ORSLEAPEPURSUTITOMIPSSIBE 
GDO + 0 
LPHDGO = LOS 

С k K ЖЕ K K X K X K K K K K K K X x 

METROD RESFX 

Ы TIVE 

€ TARGET POSITION UPDATING 

E fone G= INTGRE(THDG, (-1* AT) *PITCH*D2R) 
TGTAX=TGTA*COS (TGTHDG) 
TGTAY--TGTA*SIN(TGTHDG) 
TVELX-INTGRL(TGTVXO,TGTAX) 
TVELY-INTGRL(TGTVYO,TGTAV) 
DEN ASE C TGTX0 ,TVELX) 
VI TITORLD(TGOTY0,TVELY) 

Е TGTHDG - ATAN2(TVELX,TVELY) 

e PROP NAVY MISSILE POSITION UPDATING 
ЕВРО 1001 РЫЬО5-В} 
ВОО — INTCKE (BDO, BDDOT) 
B= INTOGRL (BO, BDOT) 
PNHDG - INTGRL(PHDG,4*BDOT) 
Е сое нр 

c PYM-INTGRL(MISSYO,VM*SIN(PNHDG)) 

E HORES UR SUIT MISSIEE 


е из 


PPHDG = ATAN2 (YT PPYM XT PEXM) 

БЕЛІ INTCRO MT а E 

ЕРҮМ INTGRL(MISSYO,VM*SIN(PPHDG 
€ 


С Е PURSUIT MISSILE 


PHDG 
ATAN2 (YT +TVELYSO. S5SATTG-LEYM, XIFTIVELA NOS A TTGZUE RN 
LPXM INTGRL (MISSXO, YM COS [LPHDG 


LPYM = INTGRL(MISSYO,VM*SIN(LPHDG 
Е k k kx х k k k k k k k k k k k k k X х 
C 
DYNAMIC 
с 
С PROP NAV MISSILE GEOMETRY UPDATE 


PNAM = 4*BDOT*PNRD 
PNR=((XT-PXM)**2 + (YT-PYM)**2)**. 

PNLOS = ATAN2(YT-PYM,XT-PXM) 
PNRD-TVELX*COS(PNLOS] *VM*COS(PNHDG-PNLOS) 


= PNLOSD = (-TVELX*SIN(PNLOS) -VM*SIN (PNHDG-PNLOS) ) /PNR 
С PURE PURSUIT GEOMETRY UPDATE 
PPR = ((XT-PPXM)**2 + (YT-PPYM)**2)**.5 
PPLOS = ATAN2(YT-PPYM,XT-PPXM) 
PPRD = TVELX*COS(PPLOS) - VM*COS(PPHDG-PPLOS) 
Е PPLOSD = -TVELX*SIN(PPLOS)/PPR 
C LEAD PURSUIT GEOMETRY UPDATE 
LPR = ((XT-LPXM)**2 + (YT- LPYM)**2)**.5 
LPLOS = ATAN2(YT-LPYM,XT-LPXM) 
LPRD = TVELX*COS(LPLOS) - VM*COS(LPHDG-LPLOS) 
LPLOSD - (-TVELX*SIN(LPLOS) - 
VM*SIN(LPHDC LPLO OII LER 
E TTG z- | ВЕХЪРКО 
C DIRECT INTERCEPT MISSILE GEOMETRY UPDATE 
RDE((YT-MISSYO)**2 + (XT-MISSXO)**2)**.5 
DXT,DYT,FON= CHECK(RD,TIME,VM,XT,YT) 
IF (PXM .GT. XT) CALL ENDRUN 
ТЕ (Ree Lee 0 O0 THEN 
NN=NN+1 
WRITE (31,50) XT,YT,PXM,PYM 
WRITE (32,50) PPXM, PPYM, LPXM, LPYM 
WRITE (33,51) TIME, PNLOS, PPLOS,LPLOS 
WRITE (34,52) TIME, PNLOSD, PPLOSD, LPLOSD 
WRITE (36,54) TIME. PNHDG, PPHDG, LPHDG 
50 FORMAT (4(F10.2, 2X) } 
51 FORMAT (F5.2,3(2X,F10.5)) 
52 FORMAT (F5.2,3(2X,F10.5)) 
54 FORMAT (F5.3,2X,3(F10.6,2X)) 
И К=10 
С IF (PNR.LT.0.1*RO) THEN 
C K-1 
C ENDIF 
ENDIF 
K-K-1 
C 
C 
SAMPLE 


С SAVE (А) 0.1,ХТ,ҮТ,РХМ, РЕМ, РРХМ, РРУМ, БРХЫ, БРУМ 
С PRINT 1l.0,XT,YT,LPXM CPYM I GDOT LEL@s IF ET p 
CONTROL , FINTIM- ó. O DB DT= 50 1 


TERMINA 
"DHDG = ATAN2(DYT-MISSYO,DXT-MISSXO) 
С СКАРН (A/A, DE=TEK618) AT 
(sc=1600,LO=0.00) YT(SC=750,LO=0.0,PO= 16000) 
GR (A/A,OV) 


C 
PXM (SC=1600, LO=0.0,AX=OMIT) , PYM(SC=750, LO=0) 
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C APH ША 
PPXM(SC=1600,LO=0.0,AX=0MIT), PPYH(SC=750, LO=0,AX=OMIT) 
Е АРН (A/A,OV) 
LPXM(SC= 1600, LO=0.0,AX= OMIT),LPYM($62750,LO-0,AX-O MIT) 

WRI (2,15) DHDG,FON 
15 FORMAT (варел SX. F5.2) 

WRITE (1,16) NN 

16 FORMAT (F4.0) 
END 
STOP 
FORTRAN 


SUEREEUTINECCHECR (RD.,TIME,VM,XT,YT,DXT,DYT,FON) 
FE ED, TIME, VM, AT, Y1T,DXT,DYT,FON,DMISS 
DMESS — VM*TIME 
ТЕ ШЕНГЕЗ5 EE RD). TREN 

ZT ZT 


DYT = YT 
FON = TIME 
ENDIE 
RETURN 
END 


TLS 


Е THIS IS A DSL PROGRAM TO FIND THE LOS AND RANGE OF A 
C TARGET IN A CONSTANT G TURN FROM A MISSILE ON A 
© DIRECTPATH 
C 
с 
INITIAL 
CONST 
G=32.2,D2R=.0175, PITCH=2.7,K2F=1.66667,I=1.0, PI=3.14159, K=0 
READ (2,10) VT,AT,THDG, „тетхо, TGTYO 
READ (2,11) MISHDG, DON 
10 FORMAT (F6.1,2X,F5.2, Ax. Ез Fiora 
11 ORMAT (F10.7,2X,F5.2) 
S MISSILE PARAMETE 
MISSXO = 0.0 
MISSYO = 0.0 
VM = 2500.0 
С TARGET PARAMETERS 


TGTHDG-THDG*D2R 
TGTV=VT*K2F 
TGTG =-AT*G 
TGTTVELXO=TGTV*SIN(TGTHDG) 
TGTTVELYO=TE TN FEOS (TGTHADE) 

€ EQUATIONS OF MOTION 

METHOD Rises 

DERIVATIVE 
TGTHDG=INTGRL (THDG,-AT*PITCH*D2K) 
TGTAX=TGTG*COS (TG THDG) 
TGTAY=-TGTG*SIN(TGTHDG) 
TV ELX=INTGRL (TGTTVELXO , TGTA®) 
TVELY=INTGRL (TETIVELTO ЖОЛЫН) 
XT=INTGRL (TGTXO,TVELX) 
YT=INTGRL(TSTYO, TVELY) 

C MISSEILECPOSTIIXODNSOPBASE 
AM = ИТКЕ МЕЕ ШЕШЕ 
YM = INTGRL(MISSYO,VM*SIN(MISHDG 


DYNAMIC 
R-(.XT-XM)**2 + (YT-YM,**2)**.5 
LOS = ATAN2(YT-YM, XT-XM) 
RD - TVELX*COS(LOS) - VM*COS(MISHDG-LOS) 
LOSI = ( TVELY*COS(LOS) - VM*SIN(MISHDG-LOS))/R 
ТЕ „К.ЕО.10) THEN 
MM = MM+1 
WRITE (37,15) XM, YM 
WRITE (38,16) TIME,LOS,LOSD,MISHDG 
1 FORLAT (2(F9.2,3X)) 
те FOR!AT (F5.2,2X,21F10.5,2%) F7 23 
ENDIF 
К-Е- 
IF (TIME .GT. DONE) CALL ENDRUN 
SAMPLE 
CONTROL FINTIM=10.0,DELT=.01 
C PRINT 1.0,XT,YT,XM,YM,LOS,R 
SAVE (D) 0.1,XT,YT,XM,YM,R,LOS 
TERMINAL 
REAT (1,26) NN 
WRITE (1,26) MM 
2€ FORMAT (F5.0) 
C GRAPE (A/D,DE=TEK618) XT(SC- 1500 ,LO= 0), YT(PO=15000) 
C GRAPE (A/D, DE=TEK618 OV) XM(SC=1500,LO=0,AX=OMIT) , YM 
С GRAPE (A/D,DE-TEK618) TIME,LOS 
Con GRAPH (A/D,DE-TEK618) TIME,R 
\ ћ 
STOP 
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APPENDIX B. THIRD ORDER SIMULATION PROGRAM LISTING 


Е PROP NAV MISSILE PROGRAM FOR THESIS 
INITIAL 


D 
E353 3) ,PR(3 
D DIMENS 


DIMENSION 
E A PLR Correo Cap 


C 
METHOD RKSFX 
CONST 2552 


A 
| 


SN 66667 


"ОФ 


* OOQOrHt.€ 
e 


< 


С 
10 FORMA Во 2, 2х, Е1О.2) 


SN +9 *»Pto0o 
жжм-- pP 


С 
S . PS 


) 
| 
) 
) 


Лл, 


САБОР ОСЛО О О СОТО OOO=0000 << Оз". 


ІІІ ODD. Or 


lH H H H H llo И ~ 

Hf H H H H H H H j O ОМ СИР СВА aan 02020) » ДУО О• • 
Eur we 
++ + 





UU. 


оо 


ІМІ 


= 


O~ 
© 


pa 
NO 


WwW 


- + 


CoCoGotobtotor.IP!HP!LE* коо Бо о 
wr 


rJ rr rJ r rr rJ J мноо но н нон н 
0009 090100 ХАРД ХА (0000000 00 (200 00 (02 
——— UJ» —— и GT —.. 
um A at er Nr e "gt gt "U 


~ ~ ~ ~ ~ ~ ~ ~ ~ 


———UJ ОЈ СУ УСОРУ оноо мном 


NRF FE 


DDD 


O 
H 
= 
H 
н 
< 
де 
OH 


o MEASUREMENT COVARIANCE MATRIX 


e) 


2 
ОУ 
RMCOV 


D 
<< 
O 
O 
< 


08 
RMCOV (1,2) 
0.0 


~ = ~ ~ 


e INITIALIZE THE BEARING MEASUREMENT NOISE COVARIANCE 
MATRIX 
SMCOV = 0.0 


Е 
IN TT ALTZATTON OF PROE NAV MISSILE CONSTANT VELOCITY, 
ZERO ACCEL 
LOS = ATAN2(TGTYO - MISSYO , TGTXO - MISSXO) 
о И Е А ТСТУО — MISSYO)**2 )**.5 
PHDGSATANZ(TGTYO4TGTVYO*TTGO-MISSYO,TGTXO-*- 
NETO *TTGO-MISSXO) 
VMXO = VM*COS(PHDG) 
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VMYO = VM*SIN(PHDG) 
в Е 
RDKP1 = -VM*COS(BPHDG-DOS) +eTGTV*SIN(THDG) Go Еее 
RDDKF1 = 0 
SKP1 = LOS 
SDKP1 = (TGTVYO/COS(LOS) - VM*SIN(PHDG-LOS))/R 
SDDKP1 = O 
BO = LOS 
ВОО = О 
С 
ea = ККР1 
КМС (2) = КОКР1 
КМС (3) = RDDKP1 
211 = SKP1 
а = SDKP1 
ë БЭ) = SDDKP1 
S 
е TARGET POSITION UPDATING 
TGTHDG = ATAN2(TVELX,TVELY) 
TGTAX = TGTA*COS(TGTHDG) 
TGTAY=(-1*TGTA] *SIN(TGTHDG) 
TVELX = INTGRL(TGTVXO,TGTAX) 
ТУЕҺҮ=ТЇМТСК1Ь(ТСТЧҮО,ТСТАҮ) 
XT = INTGRL(TGTXO,TVELX) 
YT-INTGRL(TGTYO,TVELY) 
e 
C THIRD ORDER PROP NAVPMISSIDESIOSITIONTUI I DATIN 
BDOT = INTGRL(BDO,BDDOT) 
Б = INTGRL(BO,BDOT) 
AM = U 
MVELX = INTGRL(VMXO,-AM*SIN(PNHDG)) 
MVELY = INTGRL(VMYO,AM*COS (PNHDG) ) 
РНС. = АТАНа (МУ, МЕХ, 
БА ДС M IRE 
Е PYM=INTGRLIMISS10,MVELY) 
QE: 
E THIRD ORDER PROP NAV MISSILE GEOMETRY UPDATE 


ЕМ = ((ХТ-РХМ)%%2 * (YT-PYM)**2) **.5 

LOS 2 ATANZ(YT-PYM,XT-PXM) 

RDOTM = TGTV*SIN(TGTHDG) /COS (LOS) = 
VM*COS ( PNHDG-LOS) 


KDDOTM - TGTA*COS(TGTHDG-BhOS) 


LOSD - ( TVELX*SIN(LOS) + VM 5 НСО hh 
LOSDD - TGTAY*COS(LOS)/RM 
С COMPUTE THE ERROR TERMS 
РЕЬК (1) = КМ - RKP1 
рЕЦК (21 == КЕОЛМ ЕРКЕ 


А SDEL = LOS — SKPI 
e 
CALL 
КАЕМАН (КНС, КМ, КООТМ,КОООТМ „ЕК, КОК,КОВЕК КЕРІ RDE AOO EEE 


K,DELR,PR,RMCOV,S,LOS,LOSD,LOSDD,SK,SDK,SDDK,TIME,SKP1,... 
SDKP1,SDDKP1,SDEL,SMCOV,PS,TK,RG,SG, HH) 
BDDOT=SDDKP1+10* (SDK-BDOT)+33.33333* (LOS-B) 


д U = 4 BDOT3RDOTE 
C SAVE VALUES OF THE GAIN MATRICES 
GR11=RG (1,1) 
рок е 
GR21=RG(2,1 
GR22=RG (2,2) 


TIS 


> ММО 


МРЬЕ 


On0nonnonnoututtto 
б 
> 
< 
trj 


C PRINT 


TERMINAL 


СКА 
СКАРН 
СКАРН 
СКАРН 
СКАРН 
СКАРН 
СКАРН 
СКАРН 
СКАРН 
СКАРН 
СКАРН 
СКАРН 
СКАРН 
СКАРН 
СКАРН 
СКАРН 
СКАРН 
СКАРН 
СКАРН 
СКАРН 
СКАРН 
ND 


OF 
FORTRAN 


айы Жы а-а 
о 


tz; 
2 


На СС ФС 
rr] (/) Muy 
DOM  UNP uu 
He —~ i ft itor 


| MMW tl il 


E 
ОК 
CONTROL FINTIM=10.0,DELT=.01 


S 
1 
l 
1 
222 
a 
a 
1 
i 


ОООООООФО 


С О Мн 


NDRUN 


EDU VI EUM РУМ 

| TIME, ĠR11, GRIZ2.GR2Z21,GR22.,GR31,GR32 
TIME,GS1,652,G$3 

) TIME,LOS,LOSD,LOSDD,U 


ЯК) THEN 


fom, Blase) 


O SAVE DATA FOR USE WITH GRAFAEL 
et РХМЕРҮУМ 


( 

F 

E 

E 

T 

X 

‚В, SK 

LOSD, BDOT, SDK 
SDEL 

ae Reto DOTM, KDDK 
G 
G 


, BDDOT 
si “52,653 
R11 СЕО GR21,GR22 


f 
, 
t 
е 
t 
г 
? 
f 


„КМ, ЕКР1, RK, RDOTM, RISE. RDSI RDDO TM RDDRK?1,RDDA 
1, 


Бо ав Е Sk hOSD;BDOD SDETBOSDBDAZBDDOT,SDDR 


(1, 2 MM 


FORMAT (Fé. 


SC=1600,LO= 
PH 


Y ұран, үсте, ұрайын, | о о | и Y 


ОРЕСТ] ИЕ] 20 О О О 0700 07 2> || 
НЕ S —_-—о°© 
пиши: сзи ООО Ос: 


Vi 


STATEMENTS FOR "PLOTTING WITH GRAFAEL 
КАРН 
.0),YT(SC=500,PO=16000) 


(A/A,DE-TER6198) XT 


OV) PXM (SC=1600,AX=OMIT) PYM(SC- 500) 
,DE-TEK618) TIME,LOS(SC-z.025 21) 
OV} TIME (AX=OMIT) , „В(РО=7.5, а BN rm. 1; 
TIME (AX=OMIT) SK(AX=OMIT,SC=.025,LO=-.1) 

‘DES TEK618) TIME,LOSD 

V) TIME(AX=OMIT) , BDOT(PO=7.5) 

V) TIME (AX=OMIT) , SDK (AX=OMIT) 
‘BES =TEK618) TIME, SDEL 
) TIME, RM(SC=2000 О-0. 
(AX=OMIT) ,RK(SC= раси 0,LO 
) TIME FDDOTM[SC= 500.0) 
AX=OMIT) , RDDK (AX=OMIT) 
| TIM 
) 


TIME,GR12 
TIME,GR21 
TIME,GR22 


8 
8 
8 
( 
8 
8 
8 
8 
8 
8 
8 
8 
SIT TIME BDDOT 
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SUBROUTINE KALMAN(RNG,RM,RDOTM,RDDOTM,RK,RDK,RDDK, 
RKP1,RDKP1,RDDKP1,K,DELR, PR, ВКЫСОУ, 5, БО, БОВЕ 
SKR,SDK,SDDK,TIME,SKP1,SDKP1,SDDKP1, SDEL,SMCOUS 
Р5, ТК, КО вс ИН! 
С SUBROUTINE TO ITERATE À KALMAN FILTER FOR RANGE 


C VARIABLES 
E GIVEN THE COVARIANCE MATRIX AND OBSERVATIONS 


REAL*8 RNG(3) 
RDKPl1,RDDKPÍ,DELR(3),PR(3,3],RMCOV(3,3),S(3), 
LOS LOSD ,LOSDD SK SDK, SDDK, SEP1 SDKPI BI i j ЗРЕЛИ 
SMCOV, TIME, PS(3,3),TK,TKSQ,A,B,C,D, 


* 
ж 
ж 
ж 
TEMP2 (2,2), TEMP3 (3 кн РГЕН 
ж 
* 


КМ, КОТЫ, КВОТЫ, КЕ, ГВК, КОБЕ 


TEMP1 (3), 
DÉT,SCOV,SPHIQ3 3) RG (SS sc (SO (MMS a al 
: TKSQ = TK*TK 
Co cromo IND THE NEW VALUES OF RPHI FROM THE PREVIOUS VALUES 
MATRIX 
P E E 
МЕНДЕЙ 
КРИТ (1, 3) - -BITESO 
RPHI(2,1) = :5*A*TÉSQ 
З 
RPHI(3,1) = A*TK 
ол 3) - B*TK 4 .5*A*TKSQ 
я КРИТ (3.3) - 1 * .5*B«*TKSQ 
C FIND THE PROJECTED COVARIANCE PR(K/K-1) = 
KPHI*PR(K-1/K-1)*RPHI 
- CALL ZERO(TEMP3,3) 
DO 103 He1'3 
102 N=1,3 
TEMP3(L,M)=TEMP3(L,M) + RPHI(L,N)*PR(N,M) 
103 CONTINUE 
104 CONTINUE 
C CLEAR THE OLD COVARIANCE MATRIX 
RO(PR, 
C MULTIPLY BY RPHÍ TRANSPOSE 
БОШТО МЕТ с 
о п ет = 
PR(L,M)=PR(L,M) +TEMP3 (L,N) *RPHI (M,N) 
105 CONTINUE 
106 CONTINUE 
107 CONTINUE 
с DEFINE THE Q MATRIX OF MANEUVER COVARIANCE 
QR(1,1) + 500 
Ок (1.2 - 0.0 
o 1 ЕЕЕ? 
95122120 
OR 2.3) SEN 
MS = QR(1,3) 
QR (3,2) - QR(2,3) 
IE o 
C NOW ADD TO THE COVARIANCE MATRIX 
DO 110 M-1,3 
PR(L,M) = PR(L,M) + QR(L,M) 
110 CONTINUE 
13 CONTINUE 


T20 


EE x nen PPRT TO THE COVARIANCE OF THE PREDICTED ESTIMATE 
5 FIND THE ESTIMATE OF RANGE MATRIX AT STEP K 


C ZERO A TEMPORARY MATRIX 
CALL ZERO (TEMP2,2) 
DO 121 Ь=1 


"ny 
ШОЕ РЕМ) + RMCOV(L,M) 


He 
NN 
PO 
о 
О 
=> 
H 
ЕН 
= 
= 
trj 


UE 
PESE ИЕА TENBZ2 (2,2) = TEMNP2(1,2) *TEMP2(2,1) 


TEMP2 (2,2) /DET 

т БА ВЕ 
RU "IDIEMP2C02,1)7DET 
TEMP2(1,1) /DET 

H + R) INVERSE 

Y Li 


an 
т; 
tz 
rÜ 
[ri 
СУ 
О 
С ~ ~ mm 


J i i PR(L,N) * COVhE (N M) 


RONE 3) = 0 
CONTINÚ 
RG = PH(HPH + R) (INVERSED) H 
ZERO A TEMP MATRIX FOR THE RNG MATRIX 
О O L=1 
ТЕМР1 (1) =0 
140 CONTINUE 


романа а] 

DO 141 M 
ТЕМРІ (1, 

CONTINUE 


DOT 143 N=1,3 
RNG(N) = TEMP1(N) + RNG (N) 


AAP Pe 
W WW 
ы PO 

O 
O 
= 
+3 
H 
22 
С, 
ti 


=1: 
== аа 
ПИЕТЕМР А ај + RO(L M) *  DELR (M) 


RDK = RNG (2) 
= RNG (3 
ZERO THE OLD RANGE TEMPORARY MATRIX 
DO 150 N = 1,3 
ТЕМР1 (№) = 0 
150 CONTINUE 
C FIND THE ESTIMATE OF THE STEP K+1 FOR THE RANGE MATRIX 


DO 152 L d 3 
PO 151 M P 
ТЕМР1 (1) = ТЕМР1 (1) + КРНТ(Ь,М) * ВМС(М) 
91 CONTINUE 
52 CONTINUE 


SAVE THE VALUES OF RNG(K+1/K) 
DO 153 №=1,3 


RNG(N) = TEMP1(N) 
153 CONTINUE 


ВЕЕТ I) 
КЕШКЕН 2) 
RDDKP1 NG(3) 
FIND THE COVARIANC 
I 


RG = PH(HPH+R) 
THEREFORE P(K/K) 


n t 


n m H 

я 
= 
Q 


EGO ESTIMATE 
eer GP (K/K= 1) 


СУСУ СУ 


T21 


CALL ZERO(TEMP3, 3) 
Do cd À 
Ро тел мат а 
TEMP3(L,M$ = TEMP3(L,M) + RG(L,N) *PR(N,M) 
ЖЕН CONTINUE 
MED CONTINUE 
163 CONTINUE 
L I s 
PR(L,M) = PR(L,M) - TEMP3(L,M) 
les CONTINUE 
5 С 
C NOW PR IS THE COVARIANCE OF FILTERED ESTIMATE P(K/K) 
C 
C 
C prap SUBROUTINE TO ITERATE A KALMAN FILTER FOR SIGMA 
C | GIVEN THE COVARIANCE MATRIX, OBSERVATION. 
C CALCULATE THE NEW SPHI MATRIX 
C 2 -3*RDDK/RK 
D = -3*RDK/RK 
SPLIT = 
SPHI(1, 2) = TK 
SPHI(1.3) = 5 D*TKSO 
ср, = 
d LIE = 1 4 .5*C*TKSQ 
SPHI(2, 3) = TK + .5*D*TKSO 
SPHI (3,2) = CxTK + .5*C*D*TKS 
Е SPHI(3,3) = 1 + D*TK + .5* (C+D) *TKSQ 
С CALCULATE THE NEXT PROJECTED P MATRIX FOR NEXT STEP 
C PS (K,K-1)=SPHI*PS*SPHI 
C CLEA? A TEMPORARY MATRIX 
CALL ZERO(TEMP3, 3) 
[ 3 
758 
 TEMP3(L,M)sTEMP3(L,M) * SPHI(L,N)*PS(N,M) 
2 CONTINUE 
2 CONTINUE 
Z C:NTINUE 
C ZERO THE OLD PS MATRIX 
ОНЕ АМО Ре ОН 
DO 205 L=1,3 
И: 
PS(L,M)=PS(L,M)+TEMP3(L,N) *SPHI (M,N) 
2 СБ CONTINU 
5 UE 
C DEFINE THE Q MATRIX OF MANEUVER COVARINCE 
9511,1) = 0.01 
ЕЕ 
о он = 05(1,2) 
ç (2.20) e OD 
на 
D , = ; 
62135] – 85023 
05 (3,3) = 06.01 
5 
PS(L,M) = PS(L,M) + OS(L,M) 
210 CONTINUE 


122 


0 ООООМ 


ООМ ООММ ANAAAAN adn 
NEP 


ООО 


450 


451 
С 


CONTINUE 
CALCULATE (HPH + R) INVERSE 


SCOV = PS(1,1) + SMCOV 
SG(1) = PS(1,1)/SCOV 
-- = ve de Ee 
SG(3) = PS(3,1)/SCOV 
NOW FIND THE CURRENT VALUES OF THE SIGMA MATRIX 

5(1) = S(1) + SG(1)*SDEL 
511 т Io t MEE 
S(3) = S(3) + SG(3)*SDEL 

STORE THE SIGMA MATRIX FOR USE IN THE PROGRAM 
SK = 512] 
SDK = 5(2 
SÐDK = 513) 

FIND THE NEXT VALUES OF THE SIGMA MATRIX 


РНЕ * 5 (К) 
ZERO A TEMPORARY MATRIX 


ПОРОГЕ з 
ТЕМР1(1) = 0 
CONTINUE 
Во 1 
pog И МЕТ 
TEMP1(L)z 
CONTINUE 
CONTINUE 


INPUT BACK INTO SIGMA MATRIX 


n 
Po 
TEMP1 (L) + SPHI(L,M) * S(M) 


DO 223 N= 1,3 
S(N) = TEMP1(N) 
CONTINUE 
STORE THE VALUE OF THE S MATRIX 
SKP1 = 5 (1) 
Е 
DKP1 = 
При Eco IL 3 се 3 
= = * 
II 552052 = eee ch 
PS(3,1)=PS (3.1) = Во (11) 436 (3) 
Борина SpPs(5 3) = pS(1 3| жес(2 
Р5(2,2|-Р8(2,2) Е ES(1,2)*56 (2) 
ЫС BS (2. 1. PS(1,1)*sG (9 
Е {1,3} =25 (1,3) - Р5(1,3 56 (1) 
Во = Рес (1.2) — Р5(1.2)*5С (1) 
PS(1,1)=PS(1,1) - PS(1,1)*SG(1) 
IF (HH .LE. 0.0] THEN 
HH=50 
WRITE (9,4) | TIME IS ', TIME 
WRITE (871 У ИЕ 
WRITE (91%) “KANGE PHI MATRIX' 
WRITE (9,*) (RPHI(L,M), M=1,3) 
CONTINU 
WRITE (9,*) 'SIGMA PHI MATRIX' 
DO 451 L= 


=1,3 
МЕ ЕО, БРНО М) о M=1, 3) 
CONTINUE 
SU TRE COVARTANCE MATRICES AT STEP К 
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401 


402 


ao 


WU) 
OO 
ESO 


) 'THE RANG 
(*) (PR(M,N 
) 'THE BEAR 
,*) (PS(M,N 


END 
SUBROUTINE ZERO(A,N) 
CLEAR A TEMPORARY MATRIX 


CONTINUE 
CONTINUr 
RETURN 
END 


E COVARIANCE MATRIX 
) , N=1,3) 


Je 


ING COVARIANCE MATRIX IS:' 


) , N*1,3) 
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THIRD ORDER MISSILE SIMULATION USING CONSTANT GAINS 


FOR 

C THE KALMAN FILTER. А SECOND ORDER PROP NAV REFERENCE 
€ MODEL IS SIMULATED FOR PLOTTING AND PARAMETER 
COMPARISONS. 


INITIAL 
D DIMENSION 
RNG (3) 5), „РЕЂЕ (2),05 (3), СК (3,3), ЕРНІ (3,3), 5РНІ(3,3) 
N N= 0 
C 
= 
METHOD RKSFX 
CONST G=32.2,D2R=.0175, K2F=1. 66667 
MISSXO=0.0 
МІ55Ү0-0.0 
о ое 
READ (2,10) VT,AT,THDG,TGTXO, тет YO 
10 БӘНКІШЕ(ЕСІІ 27 Е5.1 2Х,Ғ5.1,.2(2Х,Ғ10.2)) 
TGTV-VT*K2F 
TGTA--AT*G 
THDG=THDG*D2R 
TGTVXO=TGTV*SIN(THDG) 
: TGTVYO-TGTV*COS (THDG) 
С INITIALIZE, THE RANGE PHI MATRIX 
RPHI(1,2) + ТК 
RPHI(1,3) = .5*TK*TK 
RPHI(2,1) = 0 
ЫРЫТ(2,2) = 1 
RPHI(2,3) - TK 
RPHI(3,1) = 0 
RPHI (3,2) - 0 
, РАНА 3) -1 
С INITIALIZE, THE BEARING PHI MATRIX 
SPHI(1,2) = TK 
SPHI (1,3) = L5*TK*TK 
SPHI (2,1) = 0 
SP 2 2) = 1 
SPHI(2,3) = TK 
SPHI(3,1) = 0 
SPHI(3/2) = 0 
SPHI(3,3) = 1 
C CONSTANT STEADY STATE GAIN VALUES RANGE 
SETA - 0.0125 
GR (1, 3) = 0 
СЕ (21) = 0.0025 
GR (2.2) = 1.0 
GF (2 3) = 0 
GR (3,1) = 0.1250 
ЕРГЕ! 5299997 
i GR(3,3) = O 
C CONSTANT STEADY STATE GAIN VALUES BEARING 
GS (2) = 12.5 
у GS(3) = 1250.0 
C INITIALIZATION OF PROP NAV MISSILE CONSTANT VELOCITY, 
ZERO ACCEL 


LOS = ATAN2(TGTYO - MISSYO , TGTXO - MISSXO) 
Иа CS MISS X0)24 2 + (TGTYO = MISSYO)**2 )**,5 
TTGO= R/VM 
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РНрозАТАН2(10тУу0%тоту талатты E од 
X 


- vet Se eee 
VMYO = VM*SIN(PHDG 


ККР1 К 
RDKP1 = -VM*COS(PHDG-LOS) + TGTV*SIN(THDG) /COS (LOS) 


RDDKP1 = 0 
SKP1 = LOS 
SDKP1 = (TGTVYO/COS(LOS) - VM*SIN(PHDG-LOS))/R 
SDDKP1 = 0 
BO = LOS 
S BDO = 0 
ЕН = RKP1 
RNG(2) = RDKP1 
RNG(3) = RDDKP1 
s (1) = ЅКР1 
_ = SDKP1 
So = SDDKP1 
C INITIALIZATION OF SECOND ORDER PROP NAV REFERENCE 
MODEL 
GO=LOS 
Я GDO = 0 
DE V REV 
С TARGET POSITION UPDATING 
TGTHDG = ATAN2(TVELX,TVELY) 
TGTAX=TGTA*COS (TGTHDG) 
TGTAY=(-1*TGTA) *SIN (TGTHDG) 
s TNTET DT EN VISOS iO 
TVELY-INTGRL(TGTVYO,TGTAY 
XT-INTGRL(TGTXO, TVELX 
s ҮТЕ ТИСЕ (ECON EVE LY 
C 
С 


THIRD ORDER PROP NAV MISSILE POSITION UPDATING 
BDOT = INTGRL(BDO,BDDOT) 
В = INTGRL(BO, BDOT) 


AM = U 

MVELX = ДЕ E -AM*SIN(P ШІЛІГІ 
MVELY = INTGRL(VMYO,AM*COS (PNHDG) 
PNHDG = ATAN2(MVELY,MVELX) 


PXM-INTGRL(MISSXO,MVELX) 
PYM-INTGRL (MISSYO,MVELY) 


C 

E SECOND ORDER PROP NAV MISSILE 
GDDOT = -20*GDOT + (SOLOS-GAMMA)*100 
GDOT = INTGRL(GDO,GDDOT) 


GAMMA = INTGRL(GO,GDOT) 
SOHDG = INTGRL(PHDG, 4*GDOT) 
SOXM = INTGRL(MISSXÓ,VM*COS(SOHDG)) 
т SOYM = INTGRL(MISSYO,VM*SIN(SOHDG)) 
DYNAMIC 
С THIRD ORDER PROP NAV MISSILE GEOMETRY UPDATE 
RM = ((XT-PXM)**2 + (YT-PYM)**2)**.5 


LOS = АТАМ№2 (ҮТ- PYM 2T -PXM) 
RDOTM " TGTV*SIN(TGTHDG) /COS (LOS) т 
ui COSE ес 


COMPUTE THE ERROR TERMS 
БЕЛ = КМ = ВЕРУ 
DELF (2) = RDOTM - RDKP1 
SDEL = LOS —SsSHEP1 
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CALL KALMAN(RNG,RM,RDOTM,RDDOTM,RK,RDK,RDDK,RKP1, 
+RDKP1,RDDKP1,DELR,S,LOS,LOSD, LOSDD,SK,SDK, 
+SDDK, TIME, SKP1, SDKP1, SDDKP1, SDEL, GS, GR, RPHI, SPHI) 
C 


BDDOT=SDDKP1+10* (SDK-BDOT)+33.33333* (LOS-B) 
U = -4*BDOT*RDOTM 


SECOND ORDER PROP NAV MISSILE GEOMETRY UPDATE 
SOR = ((XT-SOXM)**2 + (YT-SOYM)**2) **.,5 
SOLOS = ATAN2(YT-SOYM, XT-SOXM) 

SOU = —-4*GDOT*SORD 


IF (PXM .GT. XT) THEN 
CALL ENDRUN 
ENDIF 


ao 


E PRTA FOR PLOTTING WITH DISSPLA 


T 
| ВАМ, РЕМ, БОХМ, ЗОУМ 
ШІП TIME, LOS,U 

ODR TIME SOLOS, SOU 


K=10 
IF (RW .LT. 
K=3 
ENDIF 
NN=NN+1 
ENDIF 
K=K-1 
FORMAT и 210.2)) 
FORMAT (F5.2,2X, EX 3223 114. 6) 


QOO 
~ 


1 
m 
al 

1 

1 

1 


отоо 
ососооо 


SAVE | "yos. 
) | 
) ТИ СЕТИ CRIS GR21, GR22 


B,SK 
SAVE s BDOT, SDK 
Wik, RRPL,RK,RDOTM, RDRP1,RDK, RDDOTM, RDDKP1,RDDK 


E 
STATEMENTS TO SAVE DATA FOR GRAFAEL 
SAVE SDE 
Еш КЕ ЕВВОТМ, RDDK 
И вот 
С 
ME NEDG, LOS, B,SK,LOSD,BDOT,SDK,LOSDD,BDDOT', SDDK 
ROL РАТ ЕТО. О, БВЬТ=. 01 


SAVE en А РРА, РУМ 
S1,GS2 
TERMINAL 


PRINT 
PRINT 


опопочзаппсапсопапичшчю 
а 
> 
< 
[rj 


READ (1,40) MM 
WRITE (1,40) NN 
FORMAT (F6.1) 
STATEMENTS FOR PLOTTING USING GRAFAEL 
GRAPH (A/ 


A,DE=TEK618) XT 
SC= 1600, LO=0. 200 YT (SC=500, PO=16000) 
PH (A/A,OV) PXM (SC=1600,AX=OMIT) , PYM(SC=500) 

GERPH ,DE-TEK618) TIME, Los (SC=, 025,LO=-.1) 

GRAPH 'ОУ) ТТМЕ(АХ=ОМТТ) В(РО= 7.5. SC=. 025, Lo=-: E 

Ибо ве авг 

GRAPH TM TIME Медиа. BDOT (PO= 1.8) 

АРАДА месин 

GRAPH 'РЕ=ТЕК61 TIME! НИ 2000.0, LO= 

GRAPH ,OV) TIME (AX=OMIT) ,RK(SC=2000.0 

GRAPH ,DESTEK618) TIME,RDDOTM(SC-500.0] 
GRAPH OV) TIME АХ= -OMIT),RDDK(AX-OMIT) 
H 
1 
ii 
1 


о 


` 


GRAFH , DE=TEK6 ШЕМЕ Т9 
, DE=TEK6 
, DE=TEK6 
,DE-TEKO 
,DE-TEKO6 


GRAPH 
GRAFH 
GRAPH 
GRAFH 


ГІМЕ 551 
TIME,GS2 
TIME ,GS3 
TIME OR11 


QO0OO000000000000000-—00104u 


C HiEcGQmhmtjtmuoo0otututu 
ERROR 
122 63 60 60 "1 0) 22) 20 ОООФОО > 


— йл. айлы. ЧИ жол, и sgn, yy, на. 


8 
8 
8 
( 
8 
8 
8 
8 
8 
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С GRAPH (L/H,DE-TEK618) TIME,GR12 
С GRAPH (M/H,DE-TEK618) TIME,GR21 
C GRAPH (N/H,DE=TEK618) TIME,GR22 
ENS GRAPH (O/F,DE=TEK618) TIME, BDDOT 
STOE 
FORTRAN 

SUBROUTINE 


KALMAN (RNG,RM,RDOTM,RDDOTM,RK,RDK,RDDK,RKP1,RDKP1, 
ж 


RDDKEP1,DELR,S,LOS,LOSD, LOSDD SK, DR ASDDK AMS KE IAS KE 
x SDDKP1,SDEL,GS,GR,RPHI,SPHI) 
С SUBROUTINE TO ITERATE А KALMAN FILTER FOR RANGE 
VARIABLES 
С GIVEN THE COVARIANCE MATRIX, OBSERVATIONS 
ж 


REAL*8 

КНС (3), ЕМ, ЕРОТМ, ЕРВОТМ „КК, КОК,КРОК ЕКЕ!, КОКРЕ, КОРКЕ И: 
ж 

DELR(2),8(3) ,LOS, LOSD,RPHI (3, 3) 7 Agee Dee БЕН Ы 
š GR(3,3),GS0 O) 


C FIND THE NEW VALUES OF RPHI FROM THE PREVIOUS VALUES 
OF SIGMA 
С MATRIX 


ро 110 М = 1,3 
TEMP1(N) = 0.0 
1482 CONTINUE 
C CONSTANT GAIN INPUTS 
DO 1o] 1420908 
DO 120 M=1,2 
TEMP1(L) = TEMP1(L) + GR(L,M) * DELR(M) 
ог CONTINUE 
y CONTINUE 
DO 125 N=1,3 
RNG(N) = TEMP1(N) + RNG(N) 
125 CONTINUE 
С SAVE THE VALUES OF RANGE MATRIX AT STEP K 
Ki RNG (1) 


К.К КМС (2) 
FODE RNG (3) 


ZERO THE OLD RANGE TEMPORARY MATRIX 
FIND THE ESTIMATE OF THE STEP K+1 FOR THE RANGE MATRIX 


D 5 13101 Дуко 
BC NE CU qs 

TEMP1(L) = TEMPE (GL) + ЕРНІ I Pr ph O 
CONTINUE 
CONTINUE 


SAVE THE VALUES OF RNG(K+1/K) 
DO 132 N=1 


ЕМС (№) = TEMP1(N) 
2» CONTINUE 


КЕРІ 
КОКЕ 
КӘБЕРІ 


С 
С SUBKOUTINE TO ITERATE А KALMAN FILTER FOR SIGMA 
VARIABLES 

= GIVE! THE MATRIX, OBSERVATION. 

C 

E 


өөө 
ии 


Са 
pon 


Осен 


ини 
y 
= 
Q 
Wir 


WHERE H= (1 


NOW FIND THE CURRENT VALUES OF THE SIGMA MATRIX 
SiT) - SUP US се 1 DES ERID 
5(2) - 5 (2) + GS (2) *SDEL 
$13) = 543) + CS3 DEL 
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o сасаоо 
Ne 


Обоје 


ana ан 


кю 
PO 


A THE "S Tu MATRIX FOR USE IN THE PROGRAM 
SDK = ` 
SDDK = S(3 

FIND THE NEXT VALUES OF THE SIGMA MATRIX 


S(K+1) = SPHI * S(K) 
ZERO À TEMPORARY MATRIX 


DO 140 1=1,3 
ТЕМР1(1) = 0 

CONTINUE 

DO 142 L= 


MP1(L) + SPHI(L,M) * S(M) 
CONTINUE 
CONTINUE 


INPUT BACK INTO SIGMA MATRIX 
DO 143 N= 1,3 
S(N) = TEMPÍIN) 
CONTINUE 
STORE THE VA 
SKPI 


SIRE 
SDDKP1 


| 
RETURN 
END 
SUBROUTINE ZERO(A,N) 
CLEAR A TEMPORARY MATRIX 


LUE OF THE S MATRIX 
SI 
|| 
50% 


CONTINUE 
RETURN 
END 
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